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ABSTRACT 


The  design  of  airplanes,  ships,  automobiles ,  and  so-called  ‘’sculptured 
parts"  *nvolves  the  design,  delineation,  and  mathematical  description  of 
bounding  surfaces.  A  method  is  described  which  makes  possible  the  description 
of  free-form  doubly  curved  surfaces  of  a  very'  general  kind.  An  extension  of 
these  ideas  to  hyper-surfaces  in  higher  dimensional  spaces  is  also  indicated. 

This  surface  technique  has  been  specifically  devised  for  use  in  the 
Computer-Aided  Design  Project  at  M.I.T. ,  and  has  already  been  successfully 
implemented  here  and  elsewhere. 


PRECEDING 

PAGE  BLANK 


ACKNOWLEDGEMENTS 


In  addition  to  the  invaluable  support  of  Project  MAC,  this  work  has  been 
made  possible  through  the  support  and  sponsorship  extended  to  the  M.I.T.  Elec¬ 
tronic  Systems  Laboratory  by  the  Manufacturing  Technology  Laboratory,  ASD, 
W'right-Patterson  Air  Force  Base  under  Contract  No.  AF -33(600) -42859  - 
M.I.T.  DSR  Project  8753. 

This  work  is  also  being  performed  partly  in  the  Design  Division  of  the 
M.I.T.  Mechanical  Engineering  Department  under  M.I.T.  DSR  Project  8756, 
and  partly  in  the  Aiken  Computation  Laboratory  of  Harvard  University  under  a 
contract  from  Bell  Telephone  Laboratories. 


TABLE  OF  CONTENTS 


Section  Pa; 

ABSTRACT  iii 

AKNOWLEDGEMENTS  iv 

I  INTRODUCTION  1 

II  NOTATION  7 

III  THE  SURFACE  EQUATION  9 

3.1  Boundary  Slope  Continuity  11 

3.2  Slope  Correction  Surface  13 

3.3  Higher-Order  Correction  Surfaces  15 

3.4  Matrix  Form  17 

3.5  Boundary- Curves  23 

3.6  Blending  Functions  31 

3.7  Cubic  Basis  Vector  32 

3.8  Difference  Equations  34 

IV  HYPERSURFACES  -  HIGHER  DIMENSIONS  39 

V  SURFACE  NORMAL  VECTORS  45 

5. 1  General  Slope  Continuity  Criteria  45 

5.2  Adjacent  Patch  Slope  Continuity  52 

5.3  Applications  54 

VI  CORNER  TWIST  VECTORS  59 

6.1  The  Quasi-Sphere  59 

6.2  The  Exact  Sphere  64 

Vn  RULED  SURFACES  71 

7. 1  Developable  Surfaces  71 

7.2  Plane /Surface  Intersections  74 


TABI./  OF  CONTENTS  (continued) 


Section 

vra  F  4TIONAL  POLYNOMIAL  FUNCTIONS 

Page 

75 

8.1 

Boundary  Curves  and  Blending  Functions 

75 

8.2 

Plane  Curves 

82 

8.3 

An  Example 

86 

8.4 

Plane  Curve  Through  a  Point 

89 

8.5 

Second  Derivative  Vectors 

93 

BIBLIOGRAPHY 

101 

APPENDIX 

—  Computer-Generated  Images 

103 

vi 


1 


SECTION  1 

INTRODUCTION 


The  purpose  of  this  work  is  to  present  the  mathematics  of  a  certain  class 
of  surfaces  which  are  suitable  for  the  design  and  description  of  arbitrary  shapes. 
In  tne  past,  the  subject  of  surface  mathematics  has  been  investigated,  in  ana¬ 
lytical  geometry  and  in  differential  geometry,  from  the  standpoint  of  the  analysis 
of  geometric  properties  of  surfaces  that  already  exist,  but  very  little  literature 
has  been  produced  on  the  subject  of  the  creation  of  such  surfaces.  As  a  typical 
example,  the  design  of  the  hull  of  a  racing  yacht  requires  the  description  of  a 
surface  of  considerable  subtlety'  and  complexity,  and  the  process  is  traditionally- 
carried  out  by  purely  graphical  procedures  which  are  exceedingly  laborious, 
since  they  entail  a  large  amount  of  trial  and  error  iteration  in  order  to  assure 
that  the  surface  is  completely  described,  and  is  smooth  and  "fair.”  The  design 
of  automobile  bodies  and  airplane  fuselages  is  similarly  tedious  and  time  con¬ 
suming,  although  mathematical  techniques  have  been  applied  to  aircraft  design 
for  a  number  of  years. 

A  few  papers  have  been  written  on  the  subject  of  fitting  existing  ship  hull 
shapes  by  means  of  various  types  of  polynominals,  with  the  two-fold  purpose 
of  smoothing  and  interpolating  the  information  contained  in  preliminary  graphi¬ 
cally  derived  hull  lines,  and  of  replacing  this  graphical  information  with  for¬ 
mulas  and  equations  that  will  permit  further  analytical  techniques  to  be  applied, 
such  as  structural  analysis  and  the  discipline  of  fluid  mechanics.  But  these 
mathematical  techniques  are  applicable  only  when  the  surface  has  already  been 
designed  to  some  degree  of  completeness,  so  as  to  furnish  enough  1. .formation 
to  make  the  mathematics  work. 

The  mathematical  structure  of  the  surfaces  to  be  described  in  the  follow¬ 
ing  discussion  has  been  devised  to  implement  the  surface  design  process  itself, 
so  as  to  make  it,  from  the  designer's  standpoint,  extremely  natural  and  easy. 
The  designer  himself  need  not  know  or  care  about  these  internal  mathematical 
details,  any  more  than  he  needs  to  know  the  specific  composition  of  the  pencils 
with  which  he  writes  or  the  mechanics  of  the  splines  with  which  he  now  draws 
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curves.  The  mathematics  is  relatively  simple,  but  it  is  nevertheless  too  com¬ 
plicated  for  hand  calculation,  and  is  designed  fcr  use  on  a  computer. 

In  the  design  of  a  three-dimensional  object ,  whether  it  be  an  airplane  fuse¬ 
lage,  an  automobile  body,  a  ship's  hull,  or  a  single  sculptured  part  of  a  machine, 
the  designer  requires  a  system  which  will  permit  him  to  define  a  surface  with  a 
minimum  of  input  information,  and  then  to  modify  this  surface,  if  he  feels  so 
inclined,  either  by  changing  the  original  input,  or  by  adding  more  design  con¬ 
straints  to  the  system. 

As  a  specific  example,  suppose  a  designer  wishes  to  design  an  airplane 
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fuselage,  using  the  SKETCHPAD  system.  ’  ’  He  would  like  to  be  able  to  draw 
die  outline  of  the  airplane  as  seen  from  the  side,  the  outline  of  the  airplane  as 
seen  from  above,  and  some  arbitrarily  selected  section  midships.  With  these 
three  arbitrary  curves  designed,  he  would  like  to  have  the  computer  automati¬ 
cally  and  immediately  generate  a  "fair"  surface  and  display  this  surface  to  him 
in  sufficient  detail  so  that  he  could  make  appropriate  judgments.  If  the  surface 
so  generated  does  not  satisfy  him,  he  would  perhaps  like  to  modify  his  original 
design  curves,  or  else  he  might  perhaps  like  to  add  other  new  sections  and  have 
the  computer  automatically  and  instantly  re-fair  the  surface  to  fit  this  additional 
information. 

The  following  sections  describe  a  very  simple,  flexible  and  general  class 
of  surfaces  which  are  able  to  fulfill  these  requirements.  It  will  be  shown  that 
a  single  algorithmic  structure  and  essentially  only  two  symbol  types  serve  to 
provide  the  following  features: 

1.  Smooth,  fair  surfaces  can  be  defined  by  a  minimum  number  of  curves, 
and  then  adjacent  surfaces  can  be  designed  to  match  position,  slope, 
curvature,  and  indeed  any  desired  order  of  derivative  along  the  ad¬ 
joining  boundaries. 

2.  The  design  curves  that  define  the  surface  can  be  of  any  kind  whatso¬ 
ever,  including  circles,  second-degree  curves,  polynominals,  trans¬ 
cendental,  and  also  sketched  curves  with  no  known  mathematical 
formula  whatsoever. 

3.  Some  classic  surfaces  are  not  necessarily  members  of  the  family 
of  surfaces  to  be  described;  nevertheless,  these  classic  surfaces 
can  be  matched  along  their  boundaries  to  any  order  of  derivative 
desired. 
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4.  The  arithmetic  involve  '  ;n  constructing  these  surfaces  is  extremely 
simple  and,  we  have  found,  easy  to  implement  on  a  digital  computer. 

It  also  lends  itself  to  special-purpose  computing  hardware,  such  as 
digital  or  analog  differential  analysers.  In  addition,  by  virtue  of  the 
form  of  the  algorithm,  the  parameters  that  define  the  shapes  are  ex¬ 
tremely  easy  to  compute.  (In  some  cases  they  may  require  no  com¬ 
putation  at  all. ) 

We  intend  to  develop  a  method  to  construct  complex  arbitrary  surfaces  by 
piecing  together  surface  "patches."  Each  such  patch  will  be  defined  by  four 
boundary  curves,  in  principle,  although  it  is  harmless  for  cnc  of  the  boundary 
curves  to  be  degenerate,  and  to  appear  as  a  point  instead  of  a  curve  segment. 

In  the  design  of  a  surface,  it  is  intended  that  the  designer  begin  with  a  single 
surface  patch,  or  a  very  small  number  of  patches,  and  then  subdivide  these 
regions  with  additional  design  curves  defining  boundaries  of  smaller  patches 
only  when  the  internal  surface  needs  modification.  This  is  somewhat  at  variance 
w'ith  the  customary  procedure  for  mathematical  curve  fitting  and  surface  fitting 
of  existing  curves  and  surfaces,  in  which  a  relatively  large  number  of  surface 
points  already  defined  by  some  other  procedure  are  used  to  obtain  mathematical 
expressions  for  a  surface  that  best  fits  them.  Instead,  the  system  to  be  de¬ 
scribed  is  intended  to  be  used  by  the  designer  at  the  outset,  in  the  process  of 
designing  the  surface,  ratter  than  later  on  as  a  means  for  making  it  mathe¬ 
matical. 

This  is  not  to  say  that  the  surface-patch  technique  cannot  be  used  to  for¬ 
mulate  patch-wise  mathematical  expressions  for  existing  surfaces,  but  ratter 
to  indicate  that  the  primary  purpose  of  this  surface  technique  is  to  facilitate 
the  initial  design  process  itself. 

When  the  design  process  is  completed,  the  surface  will  be  completely 
mathematically  defined,  since  this  definition  occurs  automatically  and  concur¬ 
rently  witt  design. 

Ordinarily  a  ship’s  hull  or  an  airplane  fuselage  is  described  by  certain 
important  curves  such  as,  in  the  case  of  the  hull,  a  keel  curve,  a  midships 
section,  and  a  curve  representing  the  sheer  or  deck  line;  these  curves  are 
sufficient  to  determine  a  surface,  since  they  form  the  boundaries  of  a  surface 
patch.  However,  ordinarily  this  primary  surface  will  not  have  certain  desired 
characteristics,  and  it  will  have  to  be  modified  by  introducing  additional 
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information,  such  as  for  instance  one  or  two  other  section  curves.  When  these 
additional  curves  are  introduced,  the  surface  algorithm  permits  the  computer 
to  "re-fair*'  the  original  hull  form  to  contain  these  curves. 

Similarly,  an  airplane  fuselage  can  be  designed  by  drawing  a  profile  curve, 
a  maximum  half-breadth  curve,  and  a  mid-section  of  the  fuselage.  Again  these 
curves  suffice  to  define  a  primary  surface,  which  in  most  cases  will  require 
modification  by  the  a  ddition  of  a  few  more  curves  to  make  more  explicit  the 
designer's  wishes.  A.  these  additional  curves  are  introduced,  the  original 
surface  will  be  sub-divided  into  patches,  but  the  algorithm  will  automatically 
insure  continuity  of  surface  slope  and  curvature  (if  desired)  and  will  incorporate 
these  additional  curves  into  the  surface  automatically.  This  should  make  initial 
surface  design  virtually  painless,  and  is  intended  to  remove  the  tedious  process 
of  surface  fairing  as  it  is  now  practiced  in  naval  architecture.  Airplane  fuse¬ 
lages  are  usually  somewhat  simpler  shapes  than,  say,  yacht  hulls,  and  for  a 
number  of  years  second-degree  curves  have  been  successfully  used  for  fuselage 
design;  on  the  other  hand,  naval  architects  have  steadilv  resisted  the  use  of  such 
methods  in  their  work,  since  the  complexity  of  yacht  shapes  makes  it  necessary 
to  pay  attention  to  the  irksome  details  of  the  geometry  involved,  and  second- 
degree  curves  prove  to  be  cumbersome  in  such  applications. 

The  system  that  is  described  in  this  report  is  intended  to  furnish  the  flexi¬ 
bility  that  second-degree  curve  techniques  lack,  and  to  remove  almost  entirely 
the  need  for  the  designer  to  be  an  analytical  geometer.  With  this  system  imple¬ 
mented  on  a  computer,  there  is  reason  to  believe  that  the  computer  can  take 
over  all  of  the  geometrical  and  mathematical  burden  of  the  design  process,  and 
leave  the  user  free  to  be  a  sculptor  assisted  by  an  exquisitely  skillful  mechan¬ 
ical  slave. 

Ultimately,  when  a  graphical  input-output  hardware  for  a  computer  is 
available  in  the  engineering  design  office,  these  methods  will  permit  designers 
to  delineate  complex  shapes  with  great  ease,  by  simply  drawing  the  salient 
curves  that  define  and  describe  them.  Already  experiments  along  these  lines 
are  in  progress  in  a  few  isolated  laboratories  both  in  universities  and  in 
industry'.  Very  soon  the  two  severe  handicaps  that  have  inhibited  the  wider 
use  of  such  graphical  devices  will  be  removed.  These  inhibiting  factors  have 
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been  high  cost  for  the  terminal  hardware  and  small  size  of  the  wording  area. 
Rapid  strides  are  being  made  on  both  these  fronts,  and  within  a  few  years  it  will 
be  possible  not  only  to  draw  on  a  virtually  unlimited  drawing  surface,  but  to 
draw  objects  directly  in  three-dimensional  space,  and  to  view  these  constructed 
objects  as  one  would  view  an  actual  physical  thing. 
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NOTATION 


We  shall  in  what  follows  relate  the  x,  y,  and  z  coordinates  of  points  on  a 
surface  to  two  independent  variables  u  and  w,  so  that  we  could  write 

x  =  f  (u,  w) 
y  =  g(u,  w) 
z  =  h(u,  w). 

If  the  functions  f,  g  and  h  were  specified,  then  for  a  pair  of  values  of  u  and  w, 
a  point  in  space  would  be  defined.  If  we  held  one  of  die  independent  variables 
fixed,  say  w,  then  by  allowing  u  to  vary,  the  point  in  space  would  trace  out  a 
curve.  If  subsequently  we  set  w  to  a  new  fixed  value  ami  again  allowed  u  to 
vary,  we  would  trace  out  another  curve,  and  so  on.  Clearly  by  stepping  the 
values  of  w  by  small  LrcreTuents  and  allowing  u  to  vary  after  each  such  step, 
we  could  produce  a  family  of  space  curves  that  would  lie  on  the  surface  and 
define  it.  All  that  is  needed  is  some  convenient  and  systematic  way  of  arriving 
at  the  functions  f,  g,  and  h. 

It  will  turn  out  that  the  form  of  all  of  these  three  functions  is  the  same; 
only  certain  internal  numerical  values  are  different.  In  vector  notation  we  can 
write 

[x  y  z]  *  [f<u,w)  g  (u,w)  h  (u,w)J 

Since  V  =  jx  y  zj  is  a  suitable  conventional  abbreviation  for  the  vector 
on  die  left,  we  introduce  a  similar  abbreviation  for  the  right  hand  side: 

(uw)  =  [f  (u,w)  g  (u,w)  h  (u,w)j 

Here,  in  die  abbreviated  symbol  on  the  left,  we  shall  omit  the  comma  between 
the  two  letters.  Later  on,  when  no  ambiguity  can  arise,  we  shall  omit  the 
parentheses  as  well,  and  write  simply  uw  to  stand  for  the  vector.  It  is  to  be 
remembered  that  uw  does  not  stand  for  the  ordinary  product  of  the  two  quan¬ 
tities,  but  is  merely  a  bi-literal  symbol  standing  for  a  vector  whose  components 
are  functions  of  the  two  variables. 
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We  plan  to  build  up  surfaces  by  adjoining  surface  "patches, "  in  an  analogy 
of  the  piecewise  fitt*.  of  complicated  curves  by  curve  ;  gments  suitably  joined 
together.  Accordingly,  we  shall  at  the  beginning  focus  our  attention  on  one  such 
surface  patch.  To  simplify  arithmetic,  wc  shall  stipulate  that  the  independent 
variables,  or  parameters,  u  and  w  can  take  on  only  values  between  0  and  1. 

Then  a  surface  patch  cun  be  considered  to  be  a  surface  segment  bounded  by 
four  space  curves,  (0  w),  (1  w),  (u  0)  and  (u  1). 


Here,  typically,  the  symbol  (Ow)  stands  for  the  vector  describing  the  x,  y,  and 
z  coordinates  of  points  along  the  curve  generated  by  allowing  w  to  vary  conti¬ 
nuously  from  0  to  1,  while  u  is  held  fixed  and  equal  to  0. 

We  shall  introduce  two  scalar  functions,  and  each  a  function  of  a 
single  variable.  These  will  be  referred  to  as  "blending  functions"  for  reasons 
that  will  become  clear. 

In  order  to  compress  the  surface  equation,  and  the  proofs  that  we  w  ish  to 
demonstrate,  we  shall  use  a  kind  of  indicial  notation;  we  introduce  the  indices  i 
and  j,  which  can  assume  only  the  values  P  and  1,  and  we  invoke  the  customary 
summation  convention  for  terms  with  repeated  indices.  This  convention  in  our 
case  simply  means  that  when  an  index  is  repeated  in  a  term,  we  write  out  all 
the  possible  terms  that  the  actual  indicial  values  generate,  and  then  add  them. 
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THE  SURFACE  EQUATION 

With  these  conventions  and  notational  peculiarities  in  mind,  we  write 
(uw)  -  (iw)F.(u)  *  (uj)F.(w)  -  (ij)F.(u)F.(w). 

(Typically,  the  first  term  on  the  right  expands  as  follows: 

(iw)F.(u)  =  <Ow)F0(u)  +  (lwjF^u). 

Thus  the  complete  expansion  would  consist  of  eight  terms,  if  carried  out. )  We 
shall  proceed  to  demonstrate  that  this  surface  equation  represents  a  surface 
that  contains  the  four  boundary-  curves,  and  is  thus  defined  by  them. 

We  must  make  a  stipulation,  a  weak  one,  on  the  nature  of  the  blendirg 
functions  FQ  and  F  - 

FQ(0)  =  1  Fo(l)  -  0 

F1(l)  -  1  Fx(0)  -  0 

A  further  stipulation  is  that  FQ  and  F^  be  continuous  and  monotonic  over  the 
interval. 

Now-  set  u  ~  a,  where  a  can  only  be  either  0  or  1.  Then,  substituting  in 
the  surface  equation. 

(aw)  =  (iw)F.(a)  -  (aj)F.(w)  -  (ij)F1(a)F.(w). 

Consider  F^a)  which  occurs  twice  in  the  equation.  By-  the  stipulation,  if  i  -  a, 

Fj(a)  =  1. 

Otherwise,  if  i  /  a,  F.(a)  =  0. 

Hence  all  terms  in  the  expansion  that  contain  i  ^  a  vanish;  we  can  set  i  -  a  and 
what  remains  is 

(aw)  =  (aw)F  (a)  -  (aj)F  (w)  =  (aj;F  (a)F  (w; 

a  J  ®  J 

=  (aw)  *  (aj)F.(w)  -  (a))F.(w) 

=  (aw). 
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This  stows  that  for  a  =  0  or  1,  and  hence  (aw)  =  (0w)  or  (Iw) ,  the  surface  equa¬ 
tion  reduces  to  an  identity.  This  implies  that  the  surface  contains  its  boundaries. 
An  entirely  parallel  argument  would  show  that  the  equation  also  reduces  to  an 
identity  for  the  other  two  boundaries  (uO)  and  (ul). 

Provided  a  pair  of  functions  and  are  chosen  once  anu  for  all  that 
satisfy  the  stipulations,  the  surface  equation  may  be  constructed  immediately 
and  uniquely  for  any  set  of  boundary  curves  (uO)  (ul)  (Ow)  and  (lw).  It  is  to  be 
observed  that  no  restrictions  have  been  placed  on  die  form  of  the  boundary 
curves;  there  is  perhaps  the  restriction  that  they  form  a  closed  boundary,  at 
least  at  die  corners  (ij)  =  (00),  (01),  (10),  and  (11)  otherwise  there  will  be  mul¬ 
tiple  values  within  die  surface  segment;  similarly  they  should  be  continuous 
functions,  but  apart  from  these  rather  obvious  restrictions,  they  can  be  of  any 
shape  whatever,  including  curves  that  can  only  be  represented  by  tables  of 
values. 

We  can  gain  intuitive  insight  into  the  nature  of  such  a  surface  if  we  look 
at  one  of  the  terras,  say  (iij)F.(w). 

We  have  the  expansion 

(uj)F.(w)  =  (uO)Fg(w)  +  (ui)F1(w). 

This  represents  a  weighted  average  of  die  quantities  (uO)  and  (ul).  When 
w  =  0,  F^(0j  =  1  and  F^(0)  =  0,  and  the  expression  becomes  simply  (uO).  As  w 
increases,  the  weight  of  F^(w)  decreases,  while  that  of  F^(w)  increases,  so  that 
the  surface  partakes  of  the  nature  of  both  boundary  curves.  As  w  approaches 
the  value  1,  die  influence  of  (uO)  on  the  shape  of  the  surface  gradually  disappears, 
whi’  2  die  influence  of  (ul)  gradually  becomes  dominant.  Finally,  at  w  =  1,  the 
curve  (ul)  represents  the  shape  of  the  surface.  We  can  say  that  the  surface  is 
generated  by  a  gradual  transition  from  (uO)  to  (ul),  and  that  these  two  curve 
shapes  are  ,Tblended"  together  by  virtue  of  the  blending  functions  FQ  and  F^. 

This  discussion  is  somewhat  oversimplified,  since  we  have  omitted  the  term 
(iw)F.(u)  and  it  too  plays  a  part  in  determining  the  shape  of  the  internal  surface, 
as  does  of  course  the  term  involving  the  corner  coordinates,  (ij)F.(u)F^(w). 


I 
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The*  entire  surface  equation  is  seen  to  be  symmetric  in  u  and  u .  and  In 
virtue  of  this  ami  a  secondary  symmetry  in  the  functions  F  and  F)  ,  \.t  can 
abbreviate  proofs  about  the  behavior  of  the  surface  along  all  boundaries  by 
exhibiting  a  typical  proof  for  any  one  boundary . 

3 . 1  BOUNDARY  SLOPE  CONTINUITY 

It  is  our  aim  to  design  and  delineate  complicated  surfaces  by  adjoining 
surface  patches,  in  a  piecewise  fashion.  Consider  two  such  patches  A  and  B. 


ul 


with  a  common  boundary.  For  patch  A  the  boundary'  is  (lw):  lor  patch  B  it  is 
(Ow),  and  the  vectors  of  coordinates  are  equal, 

A  (lw)  =  B  (Ow). 

Then  the  two  patches  will  be  continuous  across  their  common  boundary.  They 
will  however  in  general  be  discontinuous  in  slope  across  the  boundary ,  and  we 
wish  to  investigate  this  and  make  some  amendments  that  will  correct  this  dis¬ 
continuity  of  slope. 

We  take  the  partial  derivative  with  respect  to  u;  Our  symbolism  for  this 
partial  derivative  is  (uw)^  =  »  anc*  w*len  we  substitute,  say,  u  =  0,  we  can 

w’rite  (Ow)^  to  mean  the  value  of  the  partial  derivative  so  obtained.  Then 

(uw)^  -  (iw)  F!  (u)  +  (uj)uF.(w)  -  (ij)  F!  (u)F.(w). 

Now  substitute  u  =  a  =  0  or  1,  as  before. 

(aw)  =  (iw)  F!  (a)  +  (aj)  F  (w)  -  (ij)  F|  (a)F  (w). 

«*  A  vt  J  i  J 


■'  '"X.  .rit  ... 


I 

I 

1 


L2 
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If  we  now  place  additional  constraints  on  the  blending  functions,  that  their  first 
derivatives 

Fj(a)  =0  (a  =  either  0  cr  1) 

we  obtain  the  result 

(aw)u  =(aj)uF.(w), 
all  other  terms  vanishing. 

This  implies,  for  example,  that  when  a  =  0, 

(0w)u  *  (00)uFq(w)  +(01)uF1(w), 

or ,  the  derivative  anywhere  along  the  boundary  in  the  u  direction  (across  the 
boundary)  depends  only  upon  the  derivatives  at  the  end-points  of  the  boundary; 
it  is  entirely  independent  of  the  shapes  of  the  four  boundary  curves ,  inc’  ading 
the  boundary  (Ow)  itself. 

Thus  for  the  two  patches  A  and  B,  if 

A(10)u  •«  B(00)u 

and 

A(ll)u  =  B(01)u 

i.  e. ,  if  the  boundary  curves  are  continuous  in  slope  in  the  u  direction  at  the  ends 
of  the  contiguous  boundary  between  patches ,  we  are  guaranteed  to  have 

A(lw)u  =  B(0w)u  everywhere  along  the  boundary  regardless  of  the 
shapes  of  the  boundary  curves  of  A  and  B.  This  is  a  remarkably  power  ft  1  and 
useful  property ,  achieved  at  the  slight  expense  of  extending  the  stipulations  on 
the  Fr 

Similarly,  the  second  derivative  with  respect  to  u  is 

(uw)uu  =  <iw)  F”  (u)  +  (uj)uuF.(w)  -  (ij)  F!1  (u)F.(w) 

and  if  we  farther  stipulate  that  F!'  (a)  =  0  we  obtain 

(aw)  =  (aj)  F.(w). 

'  'uu  uu  j'  ' 
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This  establishes  second  derivative  (or  curvature)  continuity  as  an  auto¬ 
matic  and  inherent  property  of  adjacent  patches,  provided  their  boundary  curves 
have  this  kind  of  continuity  at  the  end-points  of  the  boundary.  It  is  easy  to  see 
that  we  may  escalate  in  this  way  to  any  level  of  derivative  continuity  we  wish 
along  contiguous  boundaries . 

3.2  SLOPE  CORRECTION  SURFACE 

The  surface  equation  already  described  is  very  general ,  in  the  sense  that 
it  can  contain  virtually  any  boundary  curve  we  wish,  and  it  has  certain  benign 
properties  of  derivative  matching  along  boundaries ;  nevertheless  it  is  not  a 
universal  formula  for  all  surfaces ,  and  there  are  many  that  do  not  belong  to  its 
family.  We  have  already  seen  that  surfaces  generated  by  the  surface  equation 
have  a  definite  intrinsic  slope  along  boundaries ,  whose  variation  is  rigidly 
prescribed  by  a  simple  formula.  Obviously  surfaces  exist  whose  boundary 
slopes  do  not  match  this  intrinsic  slope ,  except  at  the  end-points  of  boundaries . 
Nevertheless,  we  wish  to  be  able  to  patch  together  such  other  surfaces  with  our 
special  surfaces ,  so  as  to  have  slope  continuity  (or  continuity  of  any  level  of 
derivative). 

To  do  so,  we  introduce  a  new  surface  equation,  describing  a  slope-correc¬ 
tion  surface,  which  when  added  to  the  first  surface  equation  has  the  property  of 
leaving  the  boundaries  unchanged .  but  causing  the  derivatives  across  boundaries 
to  vary  in  any  arbitrary  way  we  wish ,  as  we  move  along  the  boundary . 

The  equation  resembles  the  first  form  very  strongly .  It  is 


(uw)  =  (iv.')uG.(u)  +(uj)wG.(w)  -  (ij)uwG.(u)G.(w). 

Here,  typically,  (iw)^  is  a  function  of  w  only,  and  describes  the  arbitrary 
variation  of  the  derivative  with  respect  to  u  as  w  varies,  along  the  curve  (iw), 
and  similarly  for  the  other  boundaries .  The  vector  (i j)  represents  the  cross 

derivatives  of  the  four  corners .  Typically , 


(00) 


uw 


dudw 


u  =  0 
w  =  0 
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The  functions  and  G.  are  again  blending  functions  or  weighting  func¬ 
tions,  but  they  have  properties  different  from  the  functions  FQ  and  F  .  We 
stipulate 

Gq(0)  ==  G^O)  =  G0(l)  -  Gjfl)  =  0 

GJ(0)  =  1  Gj(l)  =  1 

GJd)  =0  G’  (0)  =  0 

or  in  the  indicial  notation  used  earlier, 

G^a)  =  0,  a  and  i  =  0  or  1. 

GJ(a)  =  0,  a  /  i. 

G’(a)  =  1,  a  =  i. 

We  need  to  ensure  that  the  vectors  describing  the  boundaries  vanish 
identically,  and  that  the  vectors  describing  the  slope  variation  along  boundaries 
are  indeed  given  by  the  equation.  The  proof  proceeds  along  precisely  the  same 
lines  we  used  before.  First,  substitute  u  =  a.  The  equation  becomes 

(aw)  =  <iw)uG.(a)  +  (ajJ^G.fw)  -  (ij)uwG.(a)G.(w). 

=  (aj)wG,<w). 

Consider  (aj)  .  We  wish  to  have  the  correction  surface  leave  the  original 
w 

boundary  vectors  unchanged ,  and  hence  the  boundary  vectors  of  the  correction 
surface  must  vanish;  ie, 

(iw)  =  0 

(uj)  =  0 

Then  the  derivatives  of  these  boundaries  must  also  vanish;  in  particular, 

(iw)w  =  0  and  then  (ij)^  -  0,  when  w  =  j.  Hence  (aj)w  =  C. 

Thus  (aw)  =  0  indicates  the  desired  behavior  of  the  correction  surface  along  a 
boundary.  Similarly  (ua)  =  0. 
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This  demonstrates  that  the  surface  has  the  slope  variation  along  the 
boundary  as  required.  To  make  use  of  this  slope  correction  surface,  we  must 
first  determine  what  the  intrinsic  slope  of  the  surface  to  be  corrected  is ,  and 
then  we  must  subtract  this  slope  from  the  desireo  boundary  slope ,  to  yield  the 
correction  slopes  that  enter  into  the  equation.  Thus  if  (Ow)^  is  the  desired 
slope,  and 

I  <0w>u  is  the  intrinsic  slope ,  then 
C  (0w)u  will  be  the  correction  slope , 

C  (0w)u  =  (0w)u  -  I  <0w)u. 

The  correction  slopes  C  (iw)^  and  C  (uj)w  are  the  four  functions  that  enter 
into  the  slope  correction  surface.  The  desired  surface  is  obtained  by  adding 
the  correction  surface  to  the  first  surface; 

(uw)  =  I  (uw)  +  C  (uw) 

where  we  use  the  symbol  I  (uw)  to  represent  the  surface  whose  boundary  slope 
is  being  modified. 
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3.3  HIGHER-ORDER  CORRECTION  SURFACES 

Analogous  forms  may  be  obtained  for  correction  of  higher  derivatives  along 
boundaries.  For  second  derivative  correction,  the  surface  equation  is 

(uw)  •  (iw)  H.(u)  +  (uj)  H,(w)  -  (ij)  H  (u)H.(w). 

'  'm  V  '  "ww  j  "uuww  i w  j 

In  this  equation,  the  blending  functions  have  the  stipulations  that,  for  a  =  0 

or  1  as  before , 

H^a)  =  0 

Hj'(a)  =  0 

Hj"(a)  =  0,  i  ^a 
Hj'  fi)  =  1,  i  =a. 

With  these  constraints  on  the  H^,  it  is  easy  to  arrange  matters  so  that  this 
spread -order  correction  surface  is  zero  everywhere  on  the  boundary,  has  zero 
slopes  across  boundaries,  and  has  second  derivatives  across  boundaries  speci¬ 
fied  by  (iw)  and  (uj)  whatever  these  function  may  be.  The  addition  of  this 
surface  vector  to  a  given  surface  vector  will  then  provide  a  means  for  boundary 
second-derivative  correction  without  disturbing  either  the  boundary  shapes  or 
boundary  slopes. 

Although  we  have  already  carried  out  a  similar  proof  for  slope  correction, 
it  might  be  well  to  exhibit  once  again  the  course  of  the  argument. 

First,  to  show  that  the  boundary  vectors  are  zero,  substitute  u  =  a: 

(a  =  0  or  1.) 

(aw)  =(iw)  H.(a)  +  (aj)  H.(w)  -  (ij)  H.(a)H.(w) 
uu  i  '  "ww  i  '  "uuww  i  j 

-  (ai)wwHi(w»- 

The  term  (aj)  refers  to  the  second  derivative  in  the  w  sense  at  each  of  the 
'  "ww 

four  corners,  such  as,  typically,  (00).  As  in  the  case  of  slope  correction,  we 
must  have 


(iw)  =  0  along  boundaries 


THE  SURFACE  EQUATION 


17 


Then  (iw)  =  0,  (i j)  =  0,  and  in  particular  (aj)  -  0,  so  that  the 
ww  ww  ww 

equation  satisfies  the  boundary  condition. 

For  boundary  slope  vectors,  differentiate  with  respect  to  u: 

»  » 

(uw)  =  (iw)  H.  (u)  +  (Uj)  H.(w)  -  (ij)  H.  (u)H.(w) 
u  uu  i  wwu  j  '  'uuww  i  j 

Set  u  =  a: 


(aw)  =  (aj)  H.(w) . 
u  ■'  wwu  j 

V.o  wish  to  have  the  slope  vectors  vanish  along  boundaries,  so  typically 


(iw)  =  0  for  all  w. 
u 

But  then  (iw)  =  0  and  (iw)  =  0  by  taking  derivatives .  The  order  of 
uw  'uww 

differentiation  is  immaterial ,  so 


(iw)  =  (iw)  ,  and  finally  we  can  conclude  that 
uww  wwu 

(ai)wwu  =  0;  a®ain  *ke  right  and  left  hand  sides  of  the  equation  are 
in  agreement. 

Finally,  we  differentiate  again  with  respect  to  U; 

«»  it 

(uw)  =  (iw)  H.  (u)+(uj)  H.(w)  -  (ij)  H  (u)H.(w). 

uu  'uu  i  "wwuu  j'  '  'wwuu  i  '  '  j  ' 

Set  u  =  a;  only  terms  in  which  a  =  i  remain; 

(aw)  =  (aw)  H  "(a)  +  (aj)  H.(w)  -  (aj)  H  (a)H.(w) 

'  uu  uu  a  'wwuu  j  '  '  J  wwuu  a  '  '  y  ’ 

-  (aw) 

'  uu 

Again  we  have  demonstrated  an  identity .  The  escalation  to  any  level  of 
boundary  derivative  correction  vector  is  obvious . 


3.4  MATRIX  FORM 


The  surface  equation 

(uw)  =  (iw)F.(u)  +(uj)F.(w)  -  (ij)F.(u)F.(v.) 
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may  be  expanded  directly  into  matrices,  to  yield: 


(uw)  =  [uO  ul] 

V 

*  [Fo°  F,“] 

Ow 

F.w 

lw 

t 

-  [V  Fi“] 

00 

01 

r_ 

V  tu 

0 

4 

10 

1! 

FjW 

In  this  we  have  treated  the  indicial  form  term  by  term  in  a  straight¬ 
forward  way.  We  shall  in  what  follows  omit  parentheses,  since  nc  misunder¬ 
standing  can  arise.  Thus  typically  FQu  is  written  in  place  of  FQ(u)  as  a  matter 
of  convenience  and  economy.  Similarly,  typically  00  is  written  instead  o?  (00); 
the  reader  should  be  reminded  that  this  is  merely  a  compact  way  of  exhibiting 
the  x,  y,  z  coordinates  a;  point  (00). 


It  means; 


00-  jx(00),  y(00),  z(00) j  when  written  out  completely. 

The  three  vector  (matrix)  products  are  equivalent  to  the  following  three  products: 

V  Fju]  0  u0  ul 

0  0  0 

0  0  0 

4. 


r 

+  l1 

V  Fi“) 

■ 

0 

0 

0 

t 

A. 

Ow 

0 

0 

V 

lw 

0 

°J 

1F1WJ 

+  [l  F  u  F,uj  fo  0  0  fl  1 

u  *  | 

0  -00  -01  F,W 

0 

o  -10  -11  F,w 

L  j  1 

and  in  this  form  we  car.  perform  the  addition,  obtaining 

(uw)  =  jl  F()u  Fju]  0  u0  ul 

0w  -00  -01 


lw 


-10 


-11 
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It  is  slightly  more  convenient  to  rewrite  this  in  the  equivalent  form 


r  i 

■  i 

(UW)  =  -  [-1  Fqu  f^J 

0  uO  ul 

-1 

Ow  00  01 

V 

lw  10  11 

F.w 

1 

so  as  to  avoid  the  awkward  minus  signs  in  the  3x3  matrix. 


Two  facts  should  be  noted.  The  leading  row  vector  in  front  of  the  matrix 
and  the  trailing  column  vector  following  the  matrix  are  transposes  of  one 
another,  but  with  different  arguments;  the  matrix  represents  the  boundary  con¬ 
ditions  of  a  patch.  The  partition  00  01  is  redundant,  since  its  elements 

10  11 

must  agree  with  uj  and  iw  for  u  and  w  equal  to  0  or  1 . 


We  have  already  suggested  that  we  can  maintain  slope  continuity  across 
boundaries  by  suitable  stipulations  on  F  ,  and  we  have  also  already  suggested 
that  when  desired  we  can  adjust  slopes  across  boundaries  bv  a  second  additive 
vector  with  suitable  stipulations  on  its  G. .  We  shall  now  investigate  the  com¬ 
bined  form  of  the  surface  equation.  To  do  so  w'e  shall  prefix  a  symbol  to  the 
vector  uw  to  indicate  whether  we  are  talking  about  the  first  surface  equation, 
or  the  correction  surface  equation,  and  we  shall  omit  the  prefix  symbol  when 
we  are  talking  about  the  combined  form.  Thus 


uw  =  suw  +•  cuw,  with 
suw  -  the  primary  surface 
cuw  =  the  correction  surface 
uw  =  the  combination. 


Accordingly,  using  this  notational  convention,  we  will  take  derivatives,  with 
respect  to  u,  of  the  surface  equation  suw  in  order  to  determine  its  slope  vector 
in  the  u  direction. 
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suw  = 
u 


[°  V  fi’uj 


0  suO  sul 

sOw  sOO  801 

slw  slO  sll 


- 1-1  v  Fiuj  r° 


We  substitute  c  =  0,  and  obtain 


sOw  = 

[sOO 

sOlJ 

F  w 

u 

L  u 

oJ 

0 

Fiw. 

Now  consider,  for  example,  sOO^  and  the  desired  00^.  The  symbol  s00u 

refers  to  the  slope  vector  at  a  corner;  we  have  already  seen  that  at  corners  the 

correction  surface  cOO  -  0 ,  and  so  00  =  sOO  .  This  is  bourne  out  intuitively 

u  u  u  J 

fay  the  reflection  that  at  (ij)  corners,  the  two  crossing  boundary  curves  com¬ 
pletely  define  the  slopes  there;  since  this  is  so,  no  correction  of  slope  need  or 
can  be  applied. 


Hence  we  should  write 

sOw  =  ]  00 

u  L  U 


Fow 
F1W.  ' 


By  analogy  and  symmetry  we  can  write  the  remaining  three  statements  ; 

slw  =  flO  11  ]  [f„w! 

u  1  u  UJ  I  0  I 


suO  =  Tf-u  f,u]  Too 
w  lo  1  J  w 


sul  =  |F_u  F,u] 

w  l  0  1  J 
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In  order  to  obtain  a  desired  slope  vector  along  any  of  the  boundaries,  we  add 
the  correction  surface,  whose  equation  is 


r 

1 

f 

cuw  =  -  -1  G„u 

G,u 

0 

cuo 

cul 

-1 

1  0 

1  J 

w 

w 

cOw 

cOO 

cOl 

Gw 

u 

uw 

uw 

0 

clw 

ciO 

ell 

Gw 

L  u 

uw 

uw 

1 

L  _) 

As  we  have  already  remarked,  the  desired  surface  uw  is  the  sum  of  the  vectors 
suw  and  cuw.  Hence  the  correction  slope  vector,  such  as,  typically,  cuO^,  is 

cuO  =  uO  -  [f.u  F.u] 
w  w  L  0  1  J 


This  is  an  entry  in  the  correction  surface  matrix. 

Now  we  introduce  a  new  fact:  the  corner  cross  derivatives  of  the  primary 
surface  equation  are  all  zero.  To  show  this,  differentiate  the  indicial  expres¬ 
sion  first  with  respect  to  u,  then  with  respect  to  w,  and  finally  set  u  =  a,  w  =  b, 
where  a  and  b  are  as  usual  either  0  or  1 .  We  have 

uwu  =  (iw)F*  (u)  +  (uj)uFj(w)  -  (ij)F'  (u)F.(w),  and 

uwuw  =  (iw)^'  (u)'  +  (uj^Fj  <w)  -  (ij)F’  (u)F’  (w). 

Evidently  this  expression  vanishes  for  (uw)  =  (a  b) .  This  shows  that  the  corner 
"twists",  or  cross  derivatives,  of  the  original  surface  all  vanish;  it  is  a 
peculiarity  of  the  first  fundamental  surface  equation. 

Hence  we  can  assert  that 


dj  =  ij  ;  This  says  that  the  desired  twists  at  corners  are 
uw  uw 

identical  with  the  correction  surface  twists ,  since  the  fundaments  surface  has 


no  twist.  We  shall  use  this  result  to  replace  the  partition 


cOO 


uw 


clO 


uw 


cOl 


uw 


ell 


uw 


with 


00 


uw 


10 


uw 


01 


uw 


11 


uw 


l 
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We  can  rewrite  the  expression  for  cuo  as  follows: 

w 

cuo  =  |1  -F 'U  -F,u| 

w  10  1  J 

and 

cul  =  [1  -F„u  -F,u| 

w  t  0  1  J 

and  of  course,  perhaps  trivially, 

0  -  [*  -Fou  -Fl“] 

Each  of  these  matric  products  represents  an  element  of  the  top  row  of  the 
correction  surface  matrix. 

Since  the  row  matrix  Jl  -Fqu  -Fjuj  is  common  to  these  three  products, 
it  can  be  factored  out  and  introduced  into  the  matrix  j-1  Gqu  G^uj  to  yield 

[-  ['  -V  -Fi“]  V  Vi 

which  is  the  same  as  the  vector 

[-1  F0a  FjU  G0«  Gju] 

We  replace  the  elements  of  the  top  row  of  the  correction  surface  matrix 


by  the  three  matrices 

o" 
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This  causes  it  to  become  a  5  x  3  matrix,  and  we  now  iiave  the  intermediate 
result, 

cuw  =  -  [-1  F0U  FjU  G0U  Gju]  [  0  -0  !  ulw  1  [  -1 


By  similar  procedures,  we  can  write  for  the  elements  of  the  first  column 
of  the  correction  surface  matrix. 


cOw  =  fow  00  01  1  1 

u  L  u  u  uJ 

-V 

l-V 

clwu  -  K  10u  Uu]  f  1  1 


-Fow 

-FjW 


and  again  trivially,  perhaps , 

0  »  [o  0  o]  1 

-Fow 

-FjW 

When  we  factor  out  the  common  column  matrix  as  before,  and  replace  each 


entry  of  the  column  matrix 


by  the  above  expressions, 


i-s.-.  -*. «,  •  •'  Wi 
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we  obtain  the  complete  matiic  expression  for  the  correction  surface: 


-  [-1  Fou  F1°  V  °1U]  [ 0 
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If  now  wc  border  the  original  surface  equation  matrix,  it  can  be  written. 


suw  =  -  £-1  Ffu  FjU  Gqu  G^J 


0  uO  ul  00 

Ow  00  01  0  0 

lw  10  11  0  0 


0  0  0  0  0  Gqw 

0  0  0  0  0  GjW 

In  this  bordering  process,  the  value  of  the  matric  product  is  unchanged. 

Since  the  pre-  and  post-multiplicative  matrices  in  this  equation  are  the 
same  as  those  of  the  correction  surface  equation,  we  can  add  the  two  5x5 
matrices  and  pre-  and  post-multiply  by  the  two  vectors.  We  shall  perform,  in 
fact, 

uw  =  suw  +  cuw,  and  obtain 

«»—  [-i  v  v  v  °iu]  r  °  i  u°  i ui  i  u°w  i  uv ]  r  -* 
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This  is  a  general  expression  for  a  slope-matching,  slope  continuous  surface  patch 
with  entirely  arbitrary  boundaries  and  entirely  arbitrary  slopes  across  these 
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boundaries .  There  are  no  stipulations  whatever  on  the  nature  of  the  boundary 
slope  function.  The  stipulation  on  the  F  and  G  functions  have  already  been  dis¬ 
cussed. 

Now  that  we  have  constructively  arrived  at  a  general  expression  for  sur¬ 
faces  that  have  a  prescribed  boundary  vector  and  a  prescribed  boundary  slope 
vector,  it  might  be  interesting  to  apply  a  proof  to  a  conjectured  higher  order 
surface  equation  in  which  not  only  boundaries ,  boundary  slopes ,  but  also 
boundary  second  derivatives  are  vector  quantities  under  control . 

We  postulate,  therefore,  that  by  analogy  the  surface  equation  is 


uw  =  -  j^-1  FqU  FjU  Gqu  GjU  Hqu  H^uJ 
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It  represents  a  surface  patch  whose  vectors  of  coordinates ,  slope ,  and  curva¬ 
ture  as  well,  are  everywhere  arbitrary  along  its  boundaries.  The  first  column 
and  first  row  of  the  7x7  matrix  represent  these  boundary  conditions;  the  re¬ 
mainder  of  the  matrix  is  redundant,  since  the  quantities  this  partition  contains 
must  all  come  from  the  column  and  row  by  differentiation. 

We  can  test  this  equation  by  seeing  whether  it  contains  a  boundary  curve. 
To  this  end ,  set  u  =  0 ,  so  that  we  check  whether  it  contains  the  boundary  (Ow) . 
We  obtain,  invoking  the  stipulations  on  the  F,  G,  and  H  functions. 
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=  _ 
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In  the  boundary  matrix  we  have  omitted  irrelevant  terms ,  because  of  the 
zero's  in  the  pre-multiplying  vector.  We  obtain,  by  performing  the  multipli¬ 
cation. 


*-  [ 


Ow  0  0  0  0  0 


o]  M 


=  Ow,  which  is  the  hoped-for  identity. 

We  can  next  try  to  see  whether  the  equation  also  conforms  to  the  boundary 
second  derivative  conditions .  It  will  be  more  convenient  in  what  follows  to 
introduce  some  abbreviated  notation. 


[to]  =  [-1  F#u  Ftu  G#u  GjU  H0u  ly] 
and  a  similar  expression  for  [fw]. 
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Set  the  7x7  boundary  condition  matrix  equal  to  ^b|  . 

With  these  abbreviations ,  the  surface  equation  is 
uw  =  -  [fuj  ^bJ  |fwj  1 . 

We  differentiate  with  respect  to  u; 

'«■-  (H  [b1  *  M  frj)  M 


uw 


and  again; 


uw 


uu 


=-(N  H  *2  M  tBJ  *  H  KJ)[H 


We  wish  to  investigate  the  right  hand  side  of  this  equation  for  u  =0,  that  is, 

for  uw  =  Ow  .  The  blending  function  vector  and  its  derivatives  become 
uu  uu 
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0 

As  for  the  first  and  second  partial  derivatives  of  the  [b]  matrix,  all 
elements  of  j^B^J  and  |^uuj  vanish  except  for  those  in  the  top  row . 

Then 

[f”o]  [b]  =  fow  00  01  00  01  00  01  ] 

t  J  L  J  L  uu  uu  uu  uuw  UUW  UUWW  UUWW  J 

2  [fo]  [bJ  =  [  the  null  vector] 

[fo]  |b  ]  =f  0  -00  -01  -00  -01 

L  J  l  UUJ  t  UU  UU  UUW  UL 


-00  -01 
UUW  UUWW  UUWW-I 


The  sum  of  these  vectors  is  evidently 


|0w  0  0  0  0  0  o] 

*■  uu  J 


Finally, 


Ow 


uu 


=  -  [owuu  0  0  0  0  0  °]  [fw] 


=  Ow  ,  as  expected , 
uu 
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We  have  she  vr.  tha?  the  extended  surface  equation  satisfies  the  second 
derivative  boundary  conditions .  In  a  similar  way  it  can  be  shown  to  satisfy  the 
first  derivative  boundary  conditions ,  but  this  was  skipped  in  favor  of  the  proof 
lor  the  higher  derivative ,  since  the  procedure  exhibits  a  few  interesting  points . 

By  analogy  we  could  construct  matrix  products  to  represent  surfaces 
which  satisfy  even  higher  derivative  conditions  across  boundaries . 

3.5  BOUNDARY  CURVES 


It  is  often  convenient  to  use  particular  boundary  curve  functions  defined 
by  the  curve  end-points  and  end-point  tangent  vectors .  We  can  use  the  blending 
functions  themselves  to  define  such  curves.  For  example,  the  uO  boundary 
curve  can  be  described  by  the  equation 


uO 


lF,“ 


Flu 
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v] 


00 

10 

00 
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10 


u 


where  the  column  vector  contains  the  end-point  information.  We  observe  that 
the  row  vector  becomes  £l  0  0  Oj  when  u  =  0;  it  becomes  [o  1  0  oj 

when  u  =  1.  Again,  if  we  take  derivatives  of  this  row  vector  with  respect  to  u 
we  obtain  JVq’u  F^'u  C^’u  G^uj  and  this  becomes  £o  0  1  oj  for  u  =  0, 

and  it  becomes  [o  0  0  l]  for  u  =  1. 

With  this  behavior  of  the  row  vector ,  it  is  easy  to  see  that  the  equation 
does  indeed  represent  a  curve  satisfying  the  end-point  conditions. 


The  matrix  form  of  the  surface  equation  has  been  shown  to  be 


uw  =  -  j-1  Fqu 

F,u  G  u 
1  0 
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Now  when,  in  computing  uw,  we  perform  the  matrix  multiplication  f  rom 
the  left ,  we  have  for  the  second  column 


[-1  F0u  FjU  G0u  Glu] 


uO 

00 

10 

00 

u 

10 

L  u. 


=  -uO  +  FJu  Gqu  G^uj  f 


u 


-  -uO  +  uO  =  0 . 


00 
i  10 

1  00 

!  f 

i 10  I 
L  UJ 


Now  if  similarly  ul,  u0^  and  ul^  are  functions  of  the  same  kind,  their 
corresponding  column  products  vanish  just  as  in  ihe  case  of  uO.  Accordingly, 
the  resulting  product  of  the  three  matrices  has  the  form 

r 


uw 


=  -  [p  0  0  oj 


-1 


V 

FjW 

V 

LGiwJ 


=  P, 


where  P  is  the  product  of  the  row  vector  and  the  first  column  of  the  matrix,  or 

p  =  h  V  Fi“  V  Giul 


■  hu  v  v  °iu] 
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If  again  the  elements  of  the  column  vector,  Ow,  lw,  Ow  and  lw  are 

u  u 

likewise  functions  described  as  outlined,  we  can  write  typically 
Ow  =  [  00  01  00 
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Fow 
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GjW 


For  the  complete  column  vector  we  have 
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Wher.  we  substitute  this  result  for  the  column  vector,  wc  obtain  the  sur  .ace 
equation 


uw  =  [  F.u  F  u  Gnu  G  u] 


V 

FjW 


G0w 

GjW 


This  is  a  particularly  convenient  form  for  computation.  The  4x4 
matrix  contains  nothing  but  information  about  the  corner  coordinates,  corner 
slopes,  and  corner  twists;  all  entries  are  constants,  and  the  partitions  of  the 
matrix  systematically  group  these  quantities.  The  leading  row  vector  and 
the  trailing  column  vector  are  transposes  of  one  another,  (but  with  different 
arguments,  of  course. ) 

We  shall  refer  to  the  4x4  matrix  as  the  "boundary  condition"  matrix, 
and  shall  assign  to  it  the  symbol  B,  so  that  the  matric  equation  for  the  surface 
could  be  written 


r 


i 
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uw  -  f  Fflu  F1u  Gftu  G^]  B  [  FQw 

F  w 

°0W 

Yj 

It  must  be  remembered  tliat  each  of  the  entries  in  B  is  a  three-vector, 
whose  components  are  x,  y,  and  z  coordinates  and  slopes  and  twists.  This 
means  that  B  is  really  a  tensor. 

3.6  BLENDING  FUNCTIONS 

We  can  relate  the  blending  function  vector  to  a  so-called  basis  vector 
in  the  following  way.  Let  [  u  u  u  u  ]  be  a  vector  whose  elements  are  a  set 

J.  Zt  o  *x 

of  linearly  independent  functions  of  the  variable  u.  Then  we  can  postulate 
the  existence  of  a  matrix  M  such  that 

[  Fqu  Fxu  Gqu  GjUj  =  fUj  «2  u3  u,]  M. 

To  evaluate  the  M  matrix,  we  substitute  u  =  0,  u  -  1  on  both  left  and  right 
hand  sides  of  the  equation.  Then  we  take  derivatives  of  both  sides,  and 
again  substitute  u  =  0  and  u  =  1.  There  results 
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The  matrix  on  the  left  is  the  identity  matrix 


10  0  0 
0  10  0 
0  0  10 
0  0  0  1 


by  virtue 
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By  the  reasoning  of  the  last  section,  we  have  for  this  vector  basis, 


0  0  0  1 
1111 
0  0  10 
3  2  10 


from  which  we  can  obtain  the  desired  inverse 


M 


2-211 
-3  3  -2  -1 

0  0  10 

_  1  0  0  0. 


Now  we  can  write 

[FQu  FjU  GQu  GjU  ] 


u  1 1  M. 


We  shall  abbreviate  the  notation  for  the  basis  vector  in  what  follows.  We 
shall  write 

[u3  u2  u  1]  =  U 
3  2 

and  { w  w  w  1  ]  .=  W. 

The  matrice  surface  equation 

uw  =  [  Fqu  FjU  Gqu  GjU  ]  B 


now  becomes,  simply  and  compactly, 

t  t 

uw  =  U  M  B  M  W  .  {Superscript  t  means  transpose. ) 

If  u  and  W  are  cubic  basis  vectors,  then  the  surface  patch  is  the  so-called 
bi-cubic  surface.  Such  surfaces  are  very  easy  to  compute,  particularly 
since  the  basis  vector  is  so  easy  to  evaluate.  In  passing  it  is  important  to 
remark  that  the  above  compact  surface  equation  is  not  limited  to  cubics; 


Fow 

F1W 

V 

G1WJ 
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U  and  W  are  not  restricted  to  cubic  basis  vectors,  and  M  is  simply  the  matrix 
that  generates  the  appropriate  set  of  blending  functions.  Among  other  pos¬ 
sibilities,  U  and  W  might  be  higher  order  polynomial  basis  vectors;  or  they 
might  be  any  set  of  linearly  independent  functions.  Provided  the  associated 
m”1  matrix  has  an  inverse,  these  basis  vectors  are  acceptable. 

We  can  write,  for  w  held  fixed,  an  expression  for  a  u  -  varying  curve 
on  the  surface: 

uw  =  U(MBMtWt)  =  UA 

where  A  is  a  column  vector  of  constant  coefficients.  We  can  write  a  simi¬ 
lar  expression  for  u  held  fixed  and  w  varying.  The  matrix  product  MBM* 
is  the  same  in  either  case.  This  suggests  that  for  any  surface  patch  this 
product  should  be  evaluated  first;  thereafter,  we  can  either  obtain  u-varying, 
w -constant  curves  or  w-varying,  u -constant  curves  in  an  obvious  way. 

We  shall  investigate  another  basis  vector  that  is  composed  of  another 
set  of  linearly  independent  functions  (not  powers  of  u)  in  a  later  article. 

3.  8  DIFFERENCE  EQUATIONS 


If  the  basis  vectors  are  polynomial  bases,  we  can  invoke  the  techniques 
of  finite  differences  to  calculate  points  on  the  surface  patch. 


Consider  the  matrix 


L  = 


10  0  0 
110  0 
1110 
1111 
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”  a 

a 

L 

b 

= 

a  +  b 

c 

a  +  b  +  c 

.  d  . 

_  a  +  b  +  c  +  d_ 

If  a,  b,  c,  d  are  respectively  third,  second,  first  and  zero-order  differences 
3 

of  the  cubic  n  ,  then  the  column  matrix  on  the  right  of  the  equation  represents 


r 
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the  corresponding  differences  for  the  cubic  (n  +  1)  .  The  differences  for 
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(n  +  k)  are  given  by 

"  a 
b 
c 
d 

multiplications  by  the  L  matrix. 

When  n  =  0,  we  can  easily  find  that 


where  L  means  k  successive 
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b 

c 

LdJ 


6 

“6 

i 

0 


for  cubics. 


Using  this. 


k  =  [0001]  L 


6 

-6 

1 

0 


In  this  expression,  the  vector  [  0  0  0  1  ]  serves  to  select  the  bottom  element 
of  the  resulting  column  vector  after  k  multiplications  by  L. 

By  extension,  we  can  write  the  more  complete  statement 


[k3  k2  k  1]  =  [0  0  0  1]  I  k 


6  0 

-6  2 

1  -1 

to  0 
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We  shall  call  the  [0  0  0  1]  vector  1;  and  we  shall  call  the  4  x  4  matrix  N, 
so  that 

[k3  k2  k  1]  =  lLkN. 


\ 


ge}®*ttSl  .  V*f'.  ‘’“.r,  ',  •, Units  J . 
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Now  let  the  usual  parametric  variable  u  be  represented  by  u  -  kS ,  where 
k  =  0,  1,  2  ...  .  and  where  8  is  an  increment  size.  Then 


,  3  2  .3  1 2  -  -  * 

|  u  u  u  1 ]  =  [  k  k  kl] 


Call  this  last  square  matrix  A .  Then 
U  =  lLkNA. 


This  expression  states  that  we  may  step  aloiig  the  u  parameter,  in  8  incre¬ 
ments,  by  successive  multiplications  by  the  L  matrix,  and  thus  evaluate 
the  U  vector  at  these  steps.  In  order  for  u  to  go  from  0  to  1,  k  must  start 
at  0  and  go  to  -g-  ,  since  k  8  ~  1. 

We  "nn  *lso  write,  for  the  W  parameter 

W  =  1  Lk  NA 


The  surface  equation  in  difference  form  can  now  be  written  out  in 

full: 

uw  =  1  Lj  NA  M  B  M*  A1  N*  L*  l\ 

Call  the  partial  product  NA  M  B  M*  A4  N*  =  S-Q,  a  square  4x4  matrix. 
Then  LS00  =  S10  a  new  square  matrix, 

and  i 

LJ  =  S.n,  after  j  multiplications. 

Uv  JO 

We  remember  that  for  any  column  of  SQ0  the  multiplication  by  L  is  a  process 
of  cumulative  addition,  as  shown  by 


a 

a 

L 

b 

- 

a  +  b 

c 

a  +  b  +  c 

d. 

a  +  b  +  c  +  d 

We  can  write  in  general  that 

L  S..  L*  =  8.  ,  .  . 

jk  j  +  1,  k  + 1 

where  the  new  square  matrix  ia  obtained  from  the  old  by  cumulative  addition 
of  column  elements,  followed  by  cumulative  addition  of  row  elements.  These 
operations  are  furthermore  commutative,  which  means  that  we  obtain  the  same 
result  if  we  first  add  row  elements  and  then  afterward  add  column  elements: 

We  have,  finally,  that  at  u  =  0,  w  =  0 
the  surface  equation  is 

00  =  1  Sqq  i  1  and  in  general 

uw  =  IS.,  1*  where  u  =  j8,  w=kS,  and  S..  has  been  formed 

JK  JK 

from  Sqq  by  j  column  additions  and  k  row  additions. 

This  obviously  furnishes  an  extremely  simple  way  to  generate  discrete  points 
on  a  surface  patch.  The  pre-multiplier  1  has  the  effect  of  selects  the  hf-iv  'ifl 
row  of  and  similarly  the  post-multiplier  1*  has  the  effect  of  Reeling  the 

last  column  of  8 . .  The  bottom  right  hand  corner  element  of  a  is  file  value 

**  s'* 

of  the  coordinate  for  a  point  on  the  surface,  at  u  =  Jo,  ^  *-•  ko. 

CoMkter  'V  repreert.  tie  rov  ohtttoed  .Ber  J  canmlm- 
five  addition  operations  have  been  performed  02  the  columns  of  8^.  The 
right  hand  element  of  fide  vector  is  iftn  value  of  the  coordinate  at  u  «  jS, 
w  =  0.  We  can  hold  u  fixed  aid  stri>  out  successive  values  of  file  coordinate 
for  w  -  varying,  simply  by  cumulative  addition  on  this  row  vector  alone,  bn 
fills  case,  the  resulting  rigid  hand  element  is  the  marching  coordinate  value. 

T 

An  analogous  remark  can  be  made  for  the  product  1  .  This  is  a 
column  vector,  and  successive  cumulative  additions  of  its  elements  marches 
out  values  of  the  surface  coordinate  for  w  =  k8  fixed,  and  u  varying. 

Although  the  arithmetic  of  the  foregoing  difference  method  is  very 
attractive,  it  possesses  certain  drawbacks  that  must  be  made  explicit.  The 
coordinate  values  are  precise  if  aid  only  if  no  truncation  error  whatever  is 
allowed  in  the  arithmetic.  Error  is  cumulative,  and  the  least  departure  from 
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We  can  write  in  general  that 

LS,  L1  =  S.  ,  ,  , 

where  the  new  square  matrix  is  obtained  from  the  old  by  cumulative  addition 
of  column  elements,  followed  by  cumulative  addition  of  row  elements.  These 
operations  are  furthermore  commutative,  which  means  that  we  obtain  the  same 
result  if  we  first  add  row  elements  and  then  afterward  add  column  elements: 

We  have,  finally,  that  at  u  =  0,  w  =  0 
the  surface  equation  is 


00  =  1  Sqq  1  and  in  general 

uw  =  IS.,  1*  where  u  =  j8,  w  =  kS,  and  S.,  has  been  formed 
jk  1  jk 

from  by  j  column  additions  and  k  row  additions. 

This  obviously  furnishes  an  extremely  simple  way  to  generate  discrete  points 

on  a  surface  patch.  The  pre-multiplier  1  has  the  effect  of  selecting  th*  b  ’n't 

row  of  S..  and  similarly  the  post-multiplier  1*  has  the  effect  of  sftlecU ng  the 
JK 

last  column  of  S., .  The  bottom  right  hand  corner  element  of  >• .  is  the  value 
jk  *  jk 

of  the  coordinate  for  a  point  on  the  surface,  at  u  -  j  u  .  w  --  k  8 . 

Consider  1  S.„.  This  represents  the  row  vector  obtained  after  j  cumula- 

j0 

tive  addition  operations  have  been  performed  on  the  columns  of  SQ^.  T  he 
right  hand  element  of  this  vector  i<s  in'5  5  rJue  of  the  coordinate  at  u  -  jS, 
w  =  0.  We  can  hold  u  fixed  a.vi  step  out  successive  values  of  the  coordinate 
for  w  -  varying,  simply  by  cumulative  addition  on  this  row  vector  alone.  In 
this  case,  the  resulting  right  hand  ele  ment  is  the  marching  coordinate  value. 

T 

An  analogous  remark  can  be  made  for  the  product  SQk  1  .  This  is  a 
column  vector,  and  successive  cumulative  additions  of  its  elements  marches 
out  vuiues  of  the  surface  coordinate  for  w  =  k8  fixed,  and  u  varying. 

Although  the  arithmetic  of  the  foregoing  difference  method  is  very 
attractive,  it  possesses  certain  drawbacks  that  must  be  made  explicit.  The 
coordinate  values  are  precise  if  and  only  if  no  truncation  error  whatever  is 
allowed  in  the  arithmetic.  Error  is  cumulative,  and  the  least  departure  from 


Menu*  vv  -w*  - ' 


39 


SECTION  IV 


HYPERSURFACES  -  HIGHER  DIMENSIONS 


We  can  readily  extend  the  surface  equation  to  describe  hyper-surfaces 
imuHureed  in  hyper-space.  For  this  purpose  we  shall  introduce  a  slight  variant 
our  notation.  We  shall  write  typically 

u  _ 

.  =  F.u.  1  =  0  or  1. 
i  i 

This  will  be  a  standard  replacement  for  the  blending  function  notation.  The 
stipulations  on  the  F^  are  as  before,  so  that  if  u  =  a,  a  =  0  or  1,  we  can  write 


a 

i 


0  when  a  4  i 


a 

i 


1  when  a  =  i. 


For  slope  continuity  across  boundaries,  typically 


=  0  where  this  symbol  means  the  first  derivative  of  the  biending 
function  with  the  argument  =  a. 

For  higher  order  continuity  across  boundaries,  the  additional  stipulations 
on  the  blending  fucntions  are  the  same  as  for  ordinary  surfaces  with  two  degrees 
of  freedom  and  have  already  been  discussed. 


The  general  surface  equation  for  hyper-space  is,  in  indicial  form, 

v  w 


(uvw  .  .  . )  =  (ujk 


) 


j  k  * 


...  .  u  w 

(jvk  .  .  .  \  .  .  ... 
i  k 


(ijw 


U  V 


1  3 


-  (N-l)  (ijk  . 


U  V  w 

•  /.  «  1  •  •  • 
1  J  k 
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In  this  eQuatiOu,  K  is  the  uuniuei  Oi  iiidependeut  pai'auietciS  in  (u  V  w. .  ,  };'  it  is 

the  number  of  degrees  of  freedom  of  a  point  on  the  hyper-surface.  The  indices 
i,  j,  k  etc.  can  take  on  only  the  values  0  or  1. 

Let  us  proceed  to  prove  that  this  surface  contains  a  boundary,  say  for 
example  the  boundary  (u  0  0  .  .  . ).  We  hope  that  the  following  equality  holds: 


(u  0  0  .  .  . )  =  (u  j  k  .  .  .) 


0  0 

•  «  »  •  • 
J  k 


+  a  ok...)“k° 


+  <i  j  o  .  .  . ) 


u  o 

i  j* 


-  (N-l)  (i  j  k  . 


u  0  0 

i  j  k  *  ‘  * 


The  last  term  in  this  expression  is  non- vanishing  if  and  only  if  all  indices  other 
than  i  are  zero,  i.  e. ,  j  *  0,  k  =  0,  etc.  We  can  accordingly  rewrite  this  term 
as 


(N-l)  (i  0  0.  .)". 

Ne:.  consider  the  second  term  on  the  right: 

(i  0  k  .  .  • )  |  k  •  •  • 

It  is  non-vanishing  if  and  only  if  k  »  0,  etc. 

We  can  accordingly  rewrite  it  as 


(  i  0  0  .  .  .)  ” 


A  similar  consideration  applies  to 

(ijO.  .  .)U°.  .  .  which  also  becomes 
(iOO.  .  .)“. 


r*rt  ,*» 
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There  are  evidently  N-l  such  terms,  all  identical,  and  they  are  removed  by  the 
last  term.  All  that  is  left  is  the  term 


but  since  j  =  0  and  k  =  0,  we  finally  have  (u  0  0  .  .  . )  on  the  right.  This  establishes 
the  identify  >  and  the  surface  equation  has  thus  been  shown  to  contain  this  boundary 
curve.  It  is  trivial  to  show  that  the  surface  contains  all  boundaries,  and  is  defined 
by  them. 

We  can  also  show  that  the  hyper-surface  contains  boundary  surfaces 
of  lower  order.  We  shall  content  ourselves  with  the  case  for  N  =  3,  and  show 
that  it  contains  surfaces  for  N  =  2  which  are  identical  with  our  ordinary  surfaces. 
We  have 


v  w 


(u  v  w)  =  (u  j  k)  k 


4  v  k) 


u  w 
i  k 


(i  jw) 


U  V 

i  j 


-  2  (i  j  k) 


U  V  w 
ij  k 


Set  v  =  0.  Then  substituting,  and  retaining  only  non- vanishing  terms, 

v  0 

(which  means  that  j  must  be  replaced  by  0  whenever  it  occurs,  and  =  ^  =  1) 


(u  0  w)  =  (u  0  1) 


w 


4  (  i  0  w)  . 


u  w 

-2<.0k)ik 


% 


4?- 
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or 


(u  0  w)  =  <u  0  k) 


w 

k 


+  (i  0  w) 


«  *  » »  u  w 
■<IOk,ik* 

This  is  the  two-degree-of-freedom  surface 

(uw)  =  (uk)  *  +  (iw)  “  -(ik)  J1  k  . 

We  shall  next  consider  the  slope  vector  of  such  a  hyper-surface, 
partial  derivatives  with  respect  to  one  of  the  variables,  say  u,  and  get 


(uvw.  .  .)u  -  (u  jk.  .  ,)u  Jk  . 


u’  w 

+  (i  vk  .  .  .)  t  k  .  . 


+  (i  j  w  .  .  . ) 


U'  V 

i  j 


hi  this,  set  u  =  0. 


u*  V  w 

-(N-l)(ijk.  .  .  k 


(0  v  w  .  .  . )  =  (0  j  k  .  .  . )  \  ™ . 

'  u  '  u  j  k 


(i  v  k  .  .  . ) 


O’  w 
i  k 


O'  v 

+  (i  j  w  .  .  . )  .  R  . 


We  take 


O'  v  w 

-(N-l)(ijk.  .  .  k.  .  . 
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or 


V  w 

(0  V  w  .  .  .  )u  =  (0  j  k  .  .  . )  u  .  k  .  .  .  , 


since  all  other  terms  vanish  by  virtue  of  .  =0. 

l 


This  result  is  analogous  to  the  one  obtained  for  boundaries  of  ordinary 
surfaces;  it  says  that  the  slope  anywhere  on  a  boundary  is  a  function  only  of  the 
slopes  at  the  "ends”  of  the  boundary,  and  are  otherwise  independent  of  the 
boundary  shapes.  Slope  continuity  across  boundaries  is  a  consequence. 


The  hyper-surface  equation  just  developed  is  defined  by  ordinary  curves, 
or  single-degree-of  freedom  boundaries;  we  can  also  write  a  hyper-surface 
equation  for  N  degrees  of  freedom,  defined  by  boundaries  with  N-l  degrees  of 
freedom.  We  shall  exhibit  the  result  for  N  -  3: 


(u  v  w) 


=  (i  v  w) 


u 

i 


+  (u  j  w)  J 

+  (u  V  k)  * 


-  (i  j  w) 


u  V 

i  j 

u  w 


/  •  ,  .  V  w 

-<“Jk'j  k 


*«>k>Mk. 


The  proof  that  this  space  contains,  for  example,  the  boundary  subspace 
(0  v  w)  follows  the  preceding  proofs  in  principle  and  will  not  be  carried  out. 
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SURFACE  NORMAL  VECTORS 

5. 1  GENERAL  SLOPE  CONTINUITY  CRITERIA 

The  surface  normal  vector  furnishes  a  convenient  mechanism  for  the  in¬ 
vestigation  of  general  criteria  for  continuity  of  surface  slope  across  boundaries 
between  surface  patches.  It  will  be  seen  that  the  continuity  conditions  already 
established  are  much  stronger  than  are  necessary,  but  that  they  are  expedient. 

Put 


u  = 

[x 

y  z  1 

l  u 

•*u  uJ 

w  = 

[x 

y  z 

L  w 

w  w* 

for  the  tangent  vectors  of  a  surface  patch  at  some  point.  Let  us  assume  that 
another  adjacent  surface  has  a  common  boundary  curve  along  u  =  constant,  w 
varying,  so  that  W  is  common  to  both  patches.  Let  the  parameter  for  this  second 
patch  be  v,  and  for  its  tangent  vector  put 

V  =  |x  y  z  j  . 

I  V  V  vJ 
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The  two  surfaces  will  be  continuous  in  slope  across  the  boundary  at  the  point  in 
question  in  case  the  three  vectors  U,  W,  and  V  are  coplanar  there. 

The  surface  normal  vector  for  one  surface  is 


N  =  U  x  W. 

If  V  is  perpendicular  to  N, 
V  •  N  =  0. 


But  then  in  this  case,  the  three  vectors  U,  W  and  V  are  coplanar  since  they  are 
all  perpendicular  to  N. 

hi  detail,  this  gives 

N  =  UxW=  [jjjl 
t  x  y  z-i 


a  vector  whose  components  are  the  familiar  Jacobians,  and 


V  •  N  =  VN 


=  fx  y  z  ] 
L  v  v  v  J 


r  i 

_ 

j 

ss 

x  y  z 

X 

V  V  V 

J 

x  y  z 

y 

u  u  u 

j 

x  y  z 

.  z_ 

ww  w 

W  — 1 

=  0 


(This  is  the  so-called  "scalar  triple  product”  of  the  vectors. )  Thus  the  vanishing 
of  the  determinant  of  the  matrix  of  the  three  tangent  vectors  is  the  general  condi¬ 
tion  for  slope  continuity  between  two  patches,  at  any  point  on  their  common 
boundary. 


This  equation  also  shows  that  we  may  have  slope  continuity  of  surfaces 
even  though  the  curvilinear  coordinates  of  the  two  surfaces  are  not  slope  - 
continuous  across  the  boundary’. 


If  the  tangent  vectors  U  and  V  are  equal  everywhere  along  the  boundary 
curve,  the  determinant  is  sure  to  vanish;  similarly  if  the  tangent  vectors  U  and 
V  are  scalar  multiples  of  one  another,  even  when  the  scalar  multiplier  is  a 
variable  quantity. 
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U£h£KAL  uGNaimicriGN  -  lANGtNi  aUm-  AutS 

y 


4- . . . - . . . . -  x 

Suppose  that  a  surface  A  already  exists,  defined  by  the  parametric  vector 
equation 

A  =  jx(<£  /?)  y (<£  Q  )  z (<f>Q  ),] 

Let  a  be  a  curve  on  the  surface;  it  is  always  possible  to  write  the  vector  equation 
for  a  in  either  of  two  forms: 

a(4>  )  =  [x(  $  )  y(  ip  )  z(  <f>  )] 

or 

a(  G  )  =  [x(  Q )  y(  Q  )  z(  Q  )] 

Suppose  we  wish  to  attach  a  surface  B  to  surface  A,  in  such  a  way  as  to  make 
curve  a  common  to  both  surfaces,  and  suppose  furthermore  that  we  wish  to 
maintain  slope  continuity  across  this  mutual  boundary. 

We  shall  consider  curve  a  to  be  the  boundary’  (Ow)  of  the  B  surface.  We 
are  at  liberty  to  design,  arbitrarily,  a  projection  of  the  other  three  boundaries, 
(uO),  (ul),  and  (lw).  Say  for  example  that  we  design  these  curves  in  the  xy  pro¬ 
jection.  Then  the  curves  represent  the  x  and  y  components  of  their  coordinate 
vectors. 
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We  identify  the  new  independent  variable  w  with  either  <f>  or  6  along  curve 
a.  The  identification  is  arbitrary,  and  we  might  for  example  use  the  linear  form 

4>  -  4>  0 

-  =w. 

*i-*o 

With  this  identification,  <f>  and  w  and,  (by  virture  of  curve  a, )  6  are  all  related, 
so  that  knowing  any  one  of  these  quantities  enables  us  to  find  the  other  two. 

We  now  have  a  correspondence  between  curve  a  and  curve  (Ow). 

We  plan  to  construct  an  F-type  surface. 


=  [l  FqU  Fju]  0  uO  ul 


Ow  -00  -01 


lw  -10  -11 


and  then  modify  its  slope  along  (Ow)  in  order  to  match  the  A-surface  slope,  by  the 
use  of  the  G-type  slope  correction  surface 


(uw)  »  1  G  u  G,u  0  u  w  ul 

'  I  0  1  J  0  w 


Ow  -00  -01  G„w 

u  uw  uw  I  0 


lw  -10  -11 

_  U  UW’  uw 


We  shall  confine  our  attention  to  the  boundary  (Ow)  where  slope-matching  is  to 
take  place.  We  shall  elect  to  accept  the  x  and  y  components  of  the  tangent-vector 
across  (0w),  as  given  by  the  F-type  surface,  and  shall  achieve  slope  matching  by 
proper  choice  of  the  z  component  of  this  vector.  We  obtain  x  and  y  components 
of  the  tangent  vector  across  (Ow)  by  the  use  of 

(0w)u  =  (00)u  FQw  +  <01)u  FjW. 

This  is  possible  because  (00)u  and  (01)u  are  known  from  (uO)  and  (ul). 
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We  are  now  ready  to  obtain  the  missing  •;  component  of  the  tangent  vector 
across  (Ow).  We  first  compute  the  surface  normal  to  A  along  curve  a.  For  this 
purpose  we  can  use  any  one  of  the  expressions 

(<£  x  (<j>  9)q 

or  (<f>  9)$  x  a (9  )q 
°r  a (<M<£  x  ($9)q 

Each  expression  yields  a  surface  normal  N;  the  three  results  are  identical. 
We  can  evaluate  this  surface  normal  vector  at  any  point  on  (Ow)  since  we  have  a 
correspondence  between  w  and  the  variables  and  Q  . 

We  also  have  the  equation 

T 

(0w)u  N  =0.  This  is  the  familiar  condition  for  surface  slope 

continuity. 

Let  N  -  [a  b  cj  after  evaluation  at  w. 

Then  the  equation  becomes 

[x(0w)u  y(0w)u  z(0w)u]  a 

b  =0 

c. 

We  already  have  the  x  and  y  components  of  this  equation,  and  can  solve  for  the 
z  component: 


-(Ow) 


a  x(0w)  +  b  y(0w) 


This  z  component  has  a  magnitude  that  ensures  that  the  complete  vector 
<0w)u  is  coplanar  with  surface  A  at  w.  Hence  (Ow)^  is  the  desired  tangent  vector 
of  surface  .  across  (Ow). 


*«.• 
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We  next  find  the  z  tangent  vector  components  2(00)^  and  z(01)u  from  (Ow)^, 
and  use  them  in  the  equation 

z(0w)  Jj  =  z{00)uFqw  +  z(01)uF1w. 

Here  the  1  superscript  indicates  that  this  is  an  intermediate  result;  it  is  the 
intrinsic  boundary  tangent  vector  for  the  F-type  surface,  and  dees  not  yet  match 
the  z(0w)u  vector  function  obtained  from  the  A  surface. 

Accordingly,  we  must  add  to  the  F-type  surface  a  G-type  tangent  vector 
correction  surface,  so  as  to  make  the  combination  have  the  desired  slope  along 
(Ow). 

This  G-type  correction  surface  is,  as  we  have  already  shown. 


(uw)c  =  [  1  Gqu  Gju] 

0 

u0C 

w 

uic 

w 

* 

1 

0wC 

u 

-ooc 

uw 

-oic 

uw 

G0w 

r 

H* 

e*o 

-ioc 

uw 

-llc 

uw 

GjW 

The  superscript  c  indicates  that  this  is  a  correction  surface. 

Slope  correction  is  necessary  only  along  the  boundary  (Ow);  we  can  enter 
c 

the  value  for  Ow  in  the  matrix,  but  the  other  entries  must  be  looked  at  in  detail, 
u 

We  have,  for  the  slope  correction  across  (Ow), 

„  c  „  .1 

Ow  =  Ow  -  Ow  . 
u  u  u 

These  latter  two  quantities  have  already  been  found  for  the  z  component,  and  so 

0w°  is  known, 
u 

Consider  Iw  .  This  is  at  a  free  boundary,  (lw),  remote  from  (Ow),  and 
^  c  c 

we  can  set  it  equal  to  0.  Then  10  and  11  are  both  zero  also.  On  the  other 

u  u 

c  c 

hand  uO  and  ul  are  connected  to  (Ow)  at  (00)  and  (01)  and  so  we  must  specify 
w  w 

them  in  such  a  way  as  to  srHsfy  the  conditions  at  these  points.  Elsewhere,  they 


too  are  arbitrary. 
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We  write  the  G  function  expression: 
<uw)  =  [l  Gqu  Gju] 


0 

uO 

ul  1 

1  * 

w 

w 

Ow 

-00 

-01 

G.w 

u 

uw 

uw 

0 

0 

0 

0 

GjW 

— 

- 

(we  omit  the  c  super-script  temporarily. ) 

Performing  the  first  multiplication,  we  have 

(uw)  =  |0w  G  u  I  uO  -00  G„u  {  ul  -01  G„u] 

LuOl  w  uw  0  •  w  uwOJ 


V 

GjW 


Consider  the  element  uO  -  00  G  u  in  the  row  vector.  Since  uO  is 

w  uw  o  w 

a  -hitrary,  it  would  be  convenient  to  choose  it  so  as  to  make  the  entire  vector 
element  vanish.  We  therefore  write 

uO  =00  G  u. 
w  uw  0 

Then  00  =  0  as  it  should,  and 

w 


since  uO  =00  G  u, 
wu  uw  0 


00  =  00  as  it  should, 

wu  uw 


Similarly,  we  may  set 


ul  =  01  G.u. 
w  uw  0 


The  result  of  these  choices  of  uO  and  ul  is  to  reduce  the  G  equation  to 

w  w 


c  c 
uw  =  Ow  G  u. 
u  0 

This  represents  the  correction  surface  z  component  which  must  be  added  to  the 
z  component  of  the  F  surface  in  order  to  obtain  slope  continuity  along  (Ow)  between 
the  given  A  surface  and  the  designed  B  surface. 
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5.2  ADJACENT- PATCH  SLOPE  CONTINUITY 

It  is  sometimes  desirable  to  define  the  boundary  curves  for  two  adjacent 
patches  so  that  at  the  junction  between  the  curves  the  tangent  vectors  have  the 
same  direction  but  are  of  different  magnitudes.  This  is  particularly  useful  when 
the  boundary  curves  are  parametric  cubics,  because  then  the  magnitudes  of  the 
tangent-vectors  at  the  end  points  control  the  behavior  of  the  curve  segment. 


01 


As  a  specific  case,  consider  the  boundary  (0w»  common  to  two  patches;  let 
the  tangent  vectors  at  (00)  and  (01)  for  the  first  patch  be  00  and  01  ^  ,  and  let 
the  tangent  vectors  for  the  next  patch  be  00^  and  01^ 

If  the  tangent  vectors  have  the  same  direction,  they  are  scalar  multiples  of  one 


* 

4> 


another. 

or 

00 

-  m  00 

u 

and 

01 

=  n  01 

Suppose  that  the  patch  ($  w)  already  exists.  We  need  to  obtain  an  appropriate 
expression  for  (uw)  so  as  to  match  surface  slopes  across  (0w).  By  the  results 
of  the  preceeding  article,  we  can  accomplish  this  in  very  general  ways,  but  in 
our  present  case  let  us  make  a  special  requirement  on  the  tangent  vectors:  let 
us  as  suit  e  tiiat  everywhere  across  (Ow)  the  tangent  vectors  have  the  same 
direction,  and  difl*.r  only  in  magnitude.  Then  for  any  w, 

0wu  =  X  0v.  £ 


SURFACE  NORMAL  VECTORS  5' 

where-  X  is  a  scalar.  We  know  that  X  takes  the  value  m  at  *v  -0,  and  takes  the 
value  n  at  w  =1,  and  we  conclude  that  X  therefore  must  be  a  scalar  function 


With  this  relationship  between  (Ow^)  and  (0*^)  the  vectcr  cross  product  is 
always  the  null  vector: 


*)u  x  {Ow)^  =  [o  0  o] 


Hence  the  scalar  triple  product  of  (0w)y,  (0w)^,  aad  (0*)^  vanishes  for  any 
(Ow)^.  This  ensures  that  the  two  surfaces  will  be  continuous  in  slope  across 
(Ow)  for  any  shape  of  (0w>  and  for  any  X  =  X  th'.  a.  >  the  proper  behavior 
at  w  =  0  and  w  =  1. 

We  could  for  example  take 

X{w)  =  m  (1-w)  +  nw. 

This  is  a  linear  variation  of  X  with  respect  to  w.  It  has  one  disadvantage, 
however,  in  that  it  introduces  un-w anted  cross  derivatives  or  twists  at  (00) 
and  (01).  la  order  to  avoid  this,  we  might  use 

X  <w)  -  m  F0w  +  n  FjW. 

Then  the  required  siope  function  across  (Ow)  for  the  (uw)  patch  is 

0*  -  (ui  F  *  +  n  F.wl  0w  , . 

u  0  l  o 

We  can  check  to  find  the  cross  derivatix'es  introduced  by  this  relationship. 
The  cross-derivative  is  obtained  by  differentiating  with  respect  to  w,  and 


\ields 


Ow  -  <m  F'  w  +  n  F*w)  0w  . 
uw  0  1  £ 

+  tm  Fqw  +  n  FjW)  Ow^. 


At  (00),  00  =  m  00 J 

'  uw  <£w 


and  at  (01),  01  *  n  0w,  , 

ow  spw 


-  -  r-v 


1  -5*  ■w’x.wreyvs  »jw» 


A«4  U^Wl>PU.Ji!lJJ 
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This  shows  that  toe  A  function  d>es  not  introduce  additional  Mists  at  the 
corners  of  the  patch,  beyond,  of  course,  those  already  inherent  in  the  (<£\v) 
surface. 


if  the  (uw)  surface  already  exists,  defined  perhaps  by  an  F-type  equation, 
its  intrinsic  tangent  vector  across  <0w)  is  a  known  function  of  (Ow),  say  Ow  , 
where  the  1  superscript  indicates  that  it  is  an  intermediate  result.  Then,  ” 
as  before,  the  correction  of  slope  is 

0wu  *  0vu  ~  0w  *  •  The  corr*jction  surface 

is  (uw}e  =  OwCG„u. 

u  0 

When  this  correction  surface  is  added  to  the  original  surface,  the  combination 
AiU  Continuous  in  slope  with  the  (<f»)  surface  across  (Ow).  The  u  and  <*> 
curvilinear  coordinates  of  the  two  patches  w  ill  be  continuous  in  slope  across 
{0w)'  but  **»*r  *a«&ent  vectors  wall  be  different  in  magnitude. 

3  APPLICATIONS 


Let  dU  represent  a  differential  vector,  sc  that 

dU  -  [  dx  d>'  d2 j  ,  in  which 

,  d  x  d  x 

dx  *  t*"  du  +  . —  dw. 
o  U  O  w 

,  dy  d  v 

dy  -  ^  du  +  a-  dw-. 

***  “  JT  du  +  r^dw  . 

<7  U  <7\V 

If  dU  is  tangent  to  a  u- varying,  w-fixed  curve,  these  become 

dx  =  du  =  x  du 
■Jv  u 


dy  -  y  du. 
u 


dz  -  du,  since  dw  =  o. 

Hence  of  =  |x  v  z  1  du. 

1  u  *  u  u  J 
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Similarlj,  if  dV\  is  a  differential  vector  tangent  to  a  w- varying,  u-fixed 


curve, 


x  y  2  dvv. 
L  w  w  \v  J 


The  normal  differential  vector  at  a  point  of  surface  will  be  given  by 
the  vector  cross  product: 


dn  =  dU  x  dW 


v  z  z  x  x  y 
’  u  u  u  u  u  u 


•  y  z  z  x  x  y 
i  w  w  w  w  |  w  w 


du  dw. 


The  determ,  nants  that  comprise' the  elements  of  the  vector  are  the 
jabobians  J  ,  J  and  J  so  thai  we  may  write 

XV  2 

dn  =  fj  J  J  1  dudw. 

1  x  y  z  ‘ 

The  magnitude  of  cn  is  equal  to  the  differential  area  of  the  elemental  parallelo¬ 
gram  described  by  dU  and  dW.  This  magnitude  is 


dn  I  =  yj  dn  •  dn  (or  - 


M  M 


=  du  dw  v  J  2  +  J,2  +  J  2  • 

X  J  z 

From  this,  we  can  construct  an  algorithm  for  finding  surface  areas  of  patches; 
we  simply  perform  numerical  integration  of  the  expression 


•1  r  1 


-  /  / 


0  0 


./J  2  +  J  ’  +  J  2  dudw. 

v  X  V  z 


Again,  if  N  is  the  unit  normal  vector  to  the  surface  at  a  point,  then 


dn  =  N  J  dn  j ,  from  which 
N.  _  iiSL 

-  hi 

j  j  j 

_  x  y  z 

s  s'  s 


where  S  =  vJ'>+J?+J2* 
x“  V  z 
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The  quantities  ^x,  fy,  **z  are  the  coordinate  components  a,  b  and  c 
S  S  S 

of  the  surface  normal. 

If  the  surface  is  to  be  manufactured  by  milling  with  a  ball-end  cutter  of 

radius  R,  then  the  cutter-center  vector  jx  y  z  1  is  related  to  the  surface 

L  C  C  0  J 

vector  [x  y  zj  by  the  simple  expression 

['c  yc  zc!  "  ix  >’  z]  *  R  [abci- 
This  describes  a  '’parallel"  surface  spaced  a  distance  R  away  from  the 
designed  surface. 

The  normal  vector  can  also  be  used  to  calculate  volumes  enclosed  by 
surface  patches  and  planes,  as  feiiov.-s. 

Suppose  we  wish  to  calculate  the  volume  contained  between  a  surface 
patch  and  the  xz  plane. 


We  can  imagine  the  volume  broken  up  into  a  number  of  slender  prisms 
whose  axes  are  all  parallel  to  the  v  ax's  (and  perpendicular  to  the  xz  plane.) 
The  area  of  the  base  of  one  of  these  prisms  is  the  projection  of  the  small 
element  of  surface  area,  or 

dA  =  J  du  du . 

y  y 
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CORNER  TWIST  VECTORS 

6.1  THE  QUASI-SPHERE 

It  is  possible  to  choose  a  parametric  cubic  that  very  nearly  approximates 
a  circle  for  one  quadrant.  We  shall  go  into  detail  about  this  shortly;  intuition 
suggests  that  similarly  we  ought  to  be  able  to  construct  an  approximation  of  an 
octant  of  a  sphere  by  means  of  a  bicubic  surface,  bounded  by  these  approxi¬ 
mations  to  circles. 

For  the  circle  approximation,  let  us  assume  that  we  will  be  content  to 
make  the  quasi-circle  pass  through  a  point  on  the  true  circle  at  u.  1  /2.  (This 
is  not  the  best  possible  approximation,  but  it  yields  quite  good  results  and  the 
arithmetic  is  simple.) 

We  shall  assume  a  circle  of  unit 
radius,  centered  at  the  origin,  with 
end-point  values  of  the  parameter  u 
as  shown.  The  tangent  vectors  are 
symmetric,  but  have  yet  undefined 
magnitudes .. 

x 


We  have 

r  3  2 

X  =  [  u  u  u  11 
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When  u  =  1/2,  this  becomes 

x  =  7  f  1  2  4  8] 


-  ;  [‘  2  4  81 


2 

-3 

0 

1 


-2 

3 

0 

0 


-2  -  a 
3  -  2a 
a 

.  0  J 


1 

-2 

1 

0 
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"  o 

1 

a 

0 

=  -  (  -  2  +  a  +  6  -  4a  +  4a)  =  -  (4  + a) 

O  O 

y/2 

a  =  8  x  -  4.  But  at  u  =  1/2,  x  =  — -  since  it  is  a  point  on  the  circle 
(by  symmetry,  at 

4 

Hence  a  =  4  (-/2  -  1).  This  is  the  required  magnitude  of  the  two  tangent 

vectors  at  u  =  0  and  u  =  1 .  (Calculation  reveals  that  the  quasi-circle  has  a 

radius  of  about  1.00016  at  u  =  1/3  (at  ~  or  30°)  so  it  is  a  good  approximation.) 

6 

We  now  establish  a  coordinate  system  for  the  sphere,  and  show  its 
boundary  curves. 


uO  =  degenerate  curve 
boundary 
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The  boundary  curves  Ow,  Iw,  and  ui  are  all  unit  circles;  the  boundary  uO 
however  is  a  degenerate  circle,  and  appears  as  a  point. 


We  shall  first  investigate  the  z  component  of  the  uw  surface  vector 


z(uw)  =  [fqu  FjU  Gqu  Gju] 

00 

01 

00 

w 

01 

w 

10 

11 

10 

w 

11 

w 

00 

u 

01 

u 

00 

uw 

01 

uw 

10 

u 

11 

u 

10 

uw 

11 

uw 

[v  v  v  °iu]  r* 


0  0 
0  0 


vl 


0  0  0 


v  v  1  w 

0  0  0 


When  we  perform  the  first  multiplication,  we  obtain 


z(uw)  =  [fqu  +  FjU  |  0  j  0  J  -a(F0“  +  F^u)J 


But  F^u  +  FjU  =  1,  by  virtue  of  the  definition  of  the  F  functions. 


Hence 


z(uw)  =  | 


1  0  0  -a 


V 

p  ttf 


=  FQw  -  a  GJw. 
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The  curvilinear  coordinates  for  w  constant  thus  yield  constant  z;  this  implies 
that  z  is  independent  of  the  other  variable  u,  and  the  w  curves  are  plane  curves. 
They  must  of  course  be  quasi -circles. 

We  have  obtained  the  value  of  the  number  a  by  investigating  a  unit  circle. 
For  a  circle  of  radius  r,  the  tangent  vector  magnitudes  must  be  proportional  to 
r,  or  equal  to  ra.  We  can  find  these  radii  for  various  values  of  w  from  either 
the  quasi-circle  y(Ow)  or  x(lw). 


We  have 
x(lw) 


|10  11  10  11  j 


w 


and 


w. 


y(0w)  --  Foo  01  00  01  1 

t  w  wJ 


F  w 
0 


F^w 

V 

LG1W 


V 

FjW 

Cow 

GjW 


in  either  case, 

r(0w)  =  [o  1  a  0  j 


F  w 
0 

F  w 
1 

G  w 
0 

G^w 


Him  a***'4 
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For  w  fixed,  the  x  and  y  coordinates  of  a  quasi -circle  are  given  by 

*  -  [F0“  V  V  v]  r  0 


r 

ar 

l  0  J 


and 


*  •  K"  Fiu  v  Giul 


r 
0 
0 
l-ar 


where  r  is  a  function  of  w,  shown  above. 


But 


r 

ar 


0  C  0  0 


0  1  a  0 


0  a  a  0 


0  0  0  0 


F  w 
0 

FjW 

V 

giwj 


and 


r 

0 

0 

-ar 


0  0  0  0 


-a  -a  0 


V 

F  w 

G0W 

lg,w 


We  obtain  these  last  results  by  simply  writing  row’s  in  the  4x4  matrix  that 
correspond  to  elements  in  the  vector  of  the  left  side.  When  we  combine 


results,  we  have 
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x(uw)  =  ^Fqu  FjU  G^u  GjUj 


y(uw)  =  [f„u  FjU  Gqu  Gju] 


‘0 

0 

0 

0  ' 

1 

►*1 

o 

< 

0 

1 

a 

0 

0 

a 

2 

a 

0 

°0W 

0 

0 

0 

.Glw. 

"o 

1 

a 

0 

Fow' 

0 

0 

0 

0 

F.w 

0 

0 

0 

0 

Go” 

0 

-a 

2 

-a 

0 

Gw 

1,  1  . 

The  equation  for  z(uw)  has  already  been  shown. 

The  striking  thing  about  the  B  matrix  as  it  appears  in  these  equations  is 
that  is  has  non-zero  entries  in  the  bottom  right  partition. 


that 

and 


00 


uw 


01 


uw 


10 


uw 


11 


uw 


*«*>„*)  =  a 


y(10uw>  =  -*  • 


By  comparison,  we  see 


These  are  the  cross  "derivatives  at  the  corners  00  and  10.  All  other  cross 
derivatives  vrnish .  We  shall  refer  to  these  cross -derivatives  as  "twists"  of 
the  surface;  uw^  is  the  twist  vector  at  a  generalized  point  on  fee  surface. 

6.2  THE  EXACT  SPHERE 

A  bi -cubic  surface  cannot  fit  a  sphere  exactly,  and  it  would  be  interesting 
to  see  whether  by  an  appropriate  choice  of  F  and  n  functions  other  than  cubics , 
an  exact  equation  can  be  constructed. 
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v 


cos  au 


where  a 
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Then 


Again 


y 


We  can  compare  the  x  equation  with 
f  2  2 

x(u)  =  jcos  au  sin  au  G^u  GjU 

2 

-  sin  an  +  aGu 
0 

sin  au  -  sin  au  +  a  G^u,  whence 

1  2 
Gqu  =  -  (sinau  -  sin  au). 


y  =  cosau 


r  2  i 

,  2 

r 

1  cos  au 

sin  au 

1  V 

G,“] 

i 

0 

0 

-a. 

cos  au 


2 

cos 


au  -  a  GjU 


G.u  -  —  (cos  au  -  cosau). 


whence 
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We  can  easily  verify  that  these  G  functions  satisfy  the  same  stipulations 
as  the  cubic  G  functions: 

G00  -  G0l  *  0 

GjO  =  Cl  -  0 

g!o  =1  g!i  -  0 

0  0 

g|o  =  0  gJi  -  1  . 

If  we  now  use  the  same  boundary  value  matrices  as  were  used  in  the 
previous  case  of  the  quasi-sphere,  lxit  with  a  =  ~  throughout,  we  obtain  for 
the  z  component  of  the  surface  vector 
z(uw)  -  Fqw  -aGw. 

This  is,  with  the  new  F  and  G  functions 

2  2 

z(uw)  =  cos  aw  -  (cos  aw  -  cos  aw) 

=  cos  aw. 

As  before,  this  shows  that  the  z  coordinate  of  the  surface  is  independent 
of  u;  the  w  curves  are  plane  curves,  and  they  are  indeed  circles.  Their  radii 
are  given  by 


=  FjW  -  a  G  w 

2  2 
=  sin  aw  (sin  aw  -  sin  aw) 

or  r  =  sin  aw. 
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x(uw)  =  [fou  FjU  G(u  Gju] 


0 

r 

ar 

.0 


r(FjU  +  aGQu)  -sin  aw  am  au, 


and 


y(mv)  =  [Fflu  FjU  Gqu  G^] 


r 

0 

0 

-ar 


-  r(F0u  -  aGjU) 

2  2 

=  sin  aw  (cos  au  -  cos  aa  *•  cos  au) 

-  sin  aw  cos  au. 

The  resulting  parametric  equations,  when  collected,  are 
x  =  sin  aw  sin  au 
y  =  sin  aw  cos  au 
z  =  cos  aw 

and  these  are  well-known. 

This  demonstrates  that  the  sphere  is  a  special  case  of  the  general  surface 
equation,  provided  the  blending  functions  are  suitably  chosen. 

The  F  and  G  functions  are  by  inspection,  seen  to  be  linear  combinations  of 
the  linearly  independent  functions  of  u, 

£cos2au  sin6  au  cos  au  sin  auj  , 
and  this  may  be  taken  as  an  appropriate  basis  vector.  Then 

FQu  FjU  Gqu  GjuJ  =  |cos2  au  sin2  au  cosau  sinauj  J^mJ 
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6a 


where  the  M  matrix  is,  in  this  ease. 


1  0 


Incidentally,  its  inverse  is 


II 
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RULED  SURFACES 


The  locus  of  straight  lines  connecting  corresponding  points  on  two  curves 
(uO)  and  (ul)  is  a  surface.  The  lines  are  called  '’rulings"  of  the  surface. 


ul 


The  equation  for  a  generalized  line  of  the  surface  is  also  the  equation  for 
the  surface: 

(u\v)  =  (ul)w  -  (uO)w  +  (uO). 

This  is  equivalent  to 

(uw)  =  (ul)w  +  (uO)  (1-w). 

From  this  equation  we  obtain  the  derivatives 
(uw)u  *  (ul)uw  *  (u0)u  (1-w) 

(uw)w  =  (ul)  -  (uO) 

(uw)  =  (ul)  -  (uO)  . 
uw  u  u 

7.1  DEVELOPABLE  SURFACES 

A  special  case  of  such  ruled  surfaces  is  of  importance  and  interest.;  If 
the  ruled  surface  has  the  property  of  being  tangent,  along  the  rulings,  to  a 
moving  plane  which  roils  around  the  surface,  then  it  may  be  deformed  by  simple 
bending  and  flattened  out  into  a  plane.  Such  a  process  is  called  "development" 
of  the  surface.  We  call  such  surfaces  "developables"  or  "wrapped  surfaces". 

(A  sheet  of  paper  can  be  wrapped  around  the  two  curves  (uO)  and  (ul)  to  form 
the  surface.  These  surfaces  are  also  known  as  "convolutes”.) 
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The  tangent-plane  condition  can  be  established  by  showing  that  the  tangent 
vector  at  a  point  on  (u0) ,  the  tangent  vector  at  a  corresponding  point  on  (ul) , 
and  die  tangent  vectors  along  the  line  joining  these  points,  are  all  coplanar.  We 
need  to  form  the  scalar  triple  product  of  these  vectors ,  and  show  that  it 
vanishes . 


The  tangent  vectors  in  question  are  (u0)u>  (ul)u>  0*0)  w  an(*  0>»,  . 

First  observe  that,  for  a  ruled  surface, 

(uw)^  =  (ul)  -  (uO) .  This  tangent  vector  is  independent  of  w,  so 

that  (uO)  =  (ul)  .  Moreover,  the  vector  is  simply  the  line  segment  joining 
w  w 

the  two  points,  as  might  be  expected.  We  can  write,  for  the  scalar  triple 
product. 


(uO) 

u 

(ul) 

u 

(ul)  -  (uO) 


0 


where  the  notation  represents  the  determinant  of  the  matrix  of  the  three  (row) 
vectors.  If  the  determinant  vanishes  for  all  values  of  u,  die  surface  is 
developable . 


The  preceding  describes  an  analytical  test  to  ensure  that  a  ruled  surface 
is  developable.  We  shall  now  describe  a  construction  that  will  enable  us  to 
define  a  ruled  surface  by  means  of  two  space  curves.  Suppose  that  the  two 
space  curves  are  defined  by  vector  functions  of  two  different  parameters, 
u  and  £ 
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•  vl 


The  scalar  triple  product  is 

(u\ 

<£l)  -  (uO) 

If  we  consider  u  the  independent  variable,  the  equation  enables  us  to  find 
<p  for  any  value  of  u;  this  value  of  <p  determines  the  point  on  (<£l)  which  corre¬ 
sponds  to  a  point  on  (uO),  so  that  the  line  joining  these  points  is  coplanar  with 
the  tangent  to  (<f>l)  and  the  tangent  to  (uO) .  We  allow  u  to  vary,  and  obtain 
related  <f>  values;  these  values  of  4>  enable  us  to  evaluate  the  components  of  the 
vector  (<f>  1).  These  components  are  the  same  as  the  components  of  the  desired 
(ul)  vector. 

Provided  we  remember  that  the  (unspecified)  functions  of  u  and  <f>  are 
different,  the  symbolism 

(ul)  =  4>  =4,  (u) 

represents  the  statement  that  the  |x  y  z|  vector  is  the  same  for  both . 

With  the  correspondence  established  between  points  on  the  two  curves, 
we  can  write  the  equation  for  the  developable  surface, 

(uw)  =  (ul)w  +  (uO)  (1  -  w) . 

This  is  the  ruled  surface  equation,  but  with  a  special  relationship  between 
curves  (uO)  and  (ul). 

7.2  PLANE/SURFACE  INTERSECTIONS 


The  general  surface  equation  can  be  cast  in  the  form 
T 

uw  =  UBW  where  U  and  W  are  vector  functions  of  u  and  w 
respectively,  and  where  B  is  a  square  matrix  describing  the  boundary  curves. 
For  example,  we  might  be  dealing  with  the  first  F-type  surface  equation. 


in  which  these  vectors  and  the  matrix  are  explicit: 
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Since  uw  is  in  reality  a  vector  consisting  of  an  x,  ay  and  a  z  component, 

there  are  three  B  matrices  which  we  can  call  B  B  B  . 

x  y  z 

We  wish  to  find  the  intersection  of  this  surface  with  the  plane 
ax  +  by  -■  f'z  +  d  =  0. 

We  can  substitute  x  =  UB  WT,  y  =  UB  WT  and  z  =  UB  W^  into  this  equation, 

x  y  z 

and  write  the  result  in  the  form 

U  jaB  +  bB  +  cB  ]  WT  +  d  =  0. 

1  x  y  z  J 

It  is  permissible  to  interchange  the  order  of  multiplication  from  aU,  bU, 
cU  to  Ua,  Ub,  Uc  because  a  b  and  c  are  scalars.  In  this  form,  the  sum 
JaB^  +  bB^  +  cB^  j  =  S,  a  square  matrix  function  of  u  and  w,  and 
T 

USW  =  -d  is  an  equation  in  the  two  variables  u  and  w..  If  w  is 
assigned  a  fixed  value,  there  results  an  equation  in  u  which  when  solved  will 
yield  a  point  on  the  intersection  curve  of  the  surface  with  the  plane,  (provided 
of  course  such  a  point  exists  for  the  chosen  value  of  w.) 

If  the  surface  in  question  is  a  bi-cubic,  the  matrix  S  is  no  longer  a 
function  of  the  variables  u  and  w,  but  consists  of  constant  elements.  In  this 
case  the  above  procedure  reduces  co  the  solution  of  a  scries  of  cubic  equations 
in  u,  where  the  coefficients  of  the  cubics  are  determined  by  successive  fixed 
values  of  w. 

In  any  case,  if  the  spacing  of  the  w  values  is  close,  the  old  value  of  u 
just  previously  determitied  for  a  particular  choice  of  w  can  appropriately  be 
used  as  a  first  trial  solution  for  the  new  value  of  w.  Algorithms  for  the 
improvement  of  this  initial  trial  value  of  u  are  not  difficult  to  construct,  and 
will  not  be  discussed  in  detail. 

If  the  plane  is  given  by,  say,  the  equation 
x  +  d  =  0 

the  solution  procedure  is  unaltered.  Not  much  simplification  results. 


SECTION  vni 


RATIONAL  POLYNOMIAL  FUNCTIONS 

8. 1  BOUNDARY  CURVES  AND  BLENDING  FUNCTIONS 

Two  kinds  of  curves  have  for  many  years  been  traditionally  used  in  air¬ 
plane  lines  design  —  cubic  polynomials,  and  conics.  Unfortunately,  each  of 
these  curve  forms  for  itself  has  certain  drawbacks.  In  the  parametric  form, 
for  ordinary  cubics,  the  entire  shape  of  a  curve  segment  is  governed  by  end 
tangent  vectors.  Sometimes  these  end  tangent  vectors  lead  to  unwanted  hooks 
and  bulges  in  the  curve  segments.  On  the  other  hand,  conics,  although  more 
benignly  behaved,  cannot  by  their  very  nature  yield  curves  with  points  of 
inflection.  Yet  such  curves  very  often  exist  in  aircraft  shapes  —  as  for  in¬ 
stance  in  the  case  of  wing  fillets. 

Because  of  these  short-comings,  a  new  curve  type  has  been  developed. 

It  is  based  upon  rational  polynomial  functions.  It  contains  both  conics  and 
ordinary  cubics  as  special  cases,  and  provides  a  great  degree  of  generality 
and  flexibility. 

We  start  by  establishing  the  form  of  this  function. 

Let  v  be  a  vector,  so  that  for  example  v  =  fx  y  z  1J  or  v  =  fx  y  1}  or 
v  =  {x  11.  The  first  of  these  can  be  thought  of  as  the  vector  (or  matrix)  of 
coordinates  on  a  space  curve;  the  second  is  the  vector  of  coordinates  for  a 
plane  curve,  and  the  last  is  the  vector  of  a  single  coordinate.  Since  this  last 
vector  yields  the  most  genera)  case,  we  shall  begin  with  it,  and  show  how  one 
might  evaluate  a  set  of  numbers  in  a  matrix  to  define  each  of  the  parametric 
coordinates  of  a  curve. 

The  product  of  v  and  a  variable  scalar  w  is  wv  =  fwx  wj.  Here  both 
wx  and  w  are  cubic  functions  of  a  parameter,  v,  and  obviously 
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This  is  the  rstio  of  two  f*iiMr  it.  *. _  »  .  . 

~  - ®  v*wuwc  iuc  ivrm  rational  tunc- 

tloa").  We  cun  repreuent  the  two  cuhtc  polynomial.  by  the  metric  equation 


wv  s  [u3  u2  u  A' 

Since  wv  =  |wx  w] .  the  matrix  A  meet  conaiet  of  four  rowa  and  two  ^i.-- 
of  conatant  coefficients.  We  now  iroceed  to  ahow  how  theae  numbera  may  be 
found  so  as  to  define  a  coordinate  of  a  curve. 

We  shall  be  interested  in  the  end-point  coordinates  of  the  curve  at  v  -  o 
and  v  «  l.  These  coordinates  are  vQ  =  [3to  !]  and  Vj  -  ^  ij  respect¬ 
ively.  A  tangent  vector  anywhere  on  the  curve  is  clearly 


v*  =  fx'  OJ 


where  the  prime  mark  means  differentiation  with  respect  to  the  parameter  v. 
The  tangent  vectors  at  v  «  0  and  v  =  1  are  therefore  v  •  =  [x  •  0}  and 

V  "  fXl'  respectively. 

Now  differentiate  both  sides  of 


wv  =  [u3  u2  u  1]  A  and  obtain 
(wv)*  =  [3u2  2u  1  OJ  A. 

Substitution  of  u  =  0  and  u  =  l  into  these  two  expressions  yields 


w  v„ 

0  0 

0  0  0  1' 

W1V1 

1111 

<*0  V 

0  0  10 

(Wj  Vj)* 

3  2  10 

uua^-yp" 
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The  4x4  matrix  on  the  right  has  an  inverse,  and  we  may  write 


0  0  0  1 

-1 

V- 

0  0 

1111 

wi  vi 

0  0  10 

<*oV 

3  2  10 

.<*1 V. 

- 

_  _ 

2-211 

v„ 

0  0 

-3  3  -2  -1 

W  V 

1  1 

0  0  10 

w  '  v„  +  w_  V  ' 

0  0  0  0 

10  0  0 

w  *  v,  +  w,  v,' 

L  1  1  1  lj 

The  square  matrix  inverse  is  constant  and  always  the  same,  and  reappears  in 
the  algebra  so  often  that  we  shall  henceforth  call  it  the  matrix  M. 

The  matrix  equation  can  be  factored  and  r  vritten  in  the  form 


*o 

0 

0 

0 

*  «■ 

vo 

0 

wl 

0 

0 

V1 

W  * 

0 

c 

w0 

0 

V 

0 

W  * 

l 

0 

w 

lj 

.V 

The  right  hand  matrix  of  v's  represents  the  desired  end  conditions  on  the 
curve.  In  our  present  case,  it  is  of  course  a  4  x  2  matrix. 

The  middle  matrix  is  4  x  4  and  contains  the  four  numbers  [wQ  w^  w^'  w^’J. 

Any  arbitrary  set  of  four  numbers  inserted  into  this  matrix  will  serve  to  define 

a  unique  pair  of  cubic  functions  of  the  parameter  u,  from  which  x  can  be  found, 

wx 

by  using,  as  we  have  said,  the  ratio  x  =  ~ . 

Instead  of  picking  these  four  numbers  arbitrarily,  however,  we  shall 
impose  further  conditions  on  the  curve  until  enough  conditions  are  imposed 
to  define  [wQ  w^  wQ’  w  ')  uniquely. 
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Furthermore,  since  in  general 
(w v)*'  =  [(wx)" 

the  quantity  w”  is  the  second  component  of  the  vector  of  (vw)u  and  therefore  is 

associated  with  the  last  column  of  the  matrix  w  v 

0  0 

W1V1 

V  vo  +  w0  V 

V  V1  +  W1 V 

But  the  last  components  of  vQ  and  vx  are  both  1 ,  and  the  last  components  of  v  • 
and  v  1  are  both  sero.  Hence 


V  “1-6  6-1  -2)  rw0-j 


We  now  can  write,  by  combining  results: 


W0  V  “1-6  6-4  -21  p»0  (V0  -  v0) 


-2  w0-  v0' 


“l  (vi  -  V 
V  <vo  -  V  -  W0  V 
.V  <vi  -  V  *  "l  v. 


* 6  *1  (51 "  V  •  •*  "o  V  -  - "  , '  <v,  -  V  - 2  “  1 V  - 2  v  v0 


Collecting, 


w0  V  *  V-4Vl,“i<G<v,  "V2  V> 

♦  "V  <-*(*! -»o»- 
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We  now  restore  this  last  expression  to  matrix  form: 

wnvAM  ”  rw0  \  V  Wl']  Hv 


6  <vi  "  vo>  “  2  V 


“2V 

."2  -  V  . 

In  oar  present  esse,  the  matrix  on  the  right  consists  of  a  column  of  numbers 
and  a  column  of  zeros.  Hence  the  column  of  zeros  can  be  discarded,  and  the 
result  is  a  4  x  1  matrix.  On  the  left,  wQ  vQ"  is  a  scalar. 

Similarly,  we  can  find  by  analogous  algebraic  procedures  that 

vr =  *wo wi v wi]  r® - vj> + 2 v l 


*«VV 

.-2V 

Then,  writing  a  matrix  equation,  we  have,  so  far, 

[wQ  v0"  WjVj”]  =  fw0w1w0*w1*l  (P  |  Q  ] 

where  [P|Q]  represents  a  4  x  2  matrix  consisting  of  the  separate  4x1  matrices 
for  wQ  vqm  and  w^  v^",  written  side  by  side  as  columns. 

We  now  introduce  another  condition.  Let  it  be  required  that  the  curve 
pass  through  the  point  v  »  [x  1]  when  u  =  *r.  (This  value  of  u  is  of 

C  C  6 

coarse  arbitrary. ) 

This  condition  leads  to 


v  = 
c 


[1  2  4  8]  M 


wovo 


W  V 
1  1 


wo  vo  +  woV 
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By  algebraic  manipulations  similar  to  the  preceding,  we  can  rearrange  the 
equation  to  read 


8  v 

c 


This  is  an  equation  in  which  v  has  two  components,  x  and  1.  The  right  hand 

c  c 

matrix  is  a  4  x  2 .  Call  it  the  R  matrix.  Then  we  can  adjoin  these  matrices, 
to  obtain 


[»oV  wivr  8vcj  =  [w0 wi  w0'  wi,]  tpQRi* 

Now  [  P  Q  R  ]  represents  a  4  x  4  matrix;  P  and  Q  are  each  4x1  matrices,  but 
R  is  a  4  x  2  matrix  We  next  transfer  wQ  v 0"  and  v^”  to  the  right  hand 
side,  obtaining 


The  right  hand  matrix  is  now  a  4  x  4.  Provided  its  determinant  does  not  vanish, 
it  has  an  inverse,  and 


[  w0  w0’  =  (0  10  |8vj  S 

where  S  is  the  4x4  inverse  of  the  matrix 

Now  that  f  wQ  Wj  wQ'  and  w  ' )  have  been  evaluated,  the  curve  is  com- 

wv 

pletely  defined,  since  the  rational  function  -  is  completely  defined. 
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8.2  PLANE  CURVES 

With  some  loss  of  generality  and  flexibility,  we  can  have  v  =  { x  y  1] , 
a  plane  curve.  We  shall  next  show  that  for  an  appropriate  choice  of 
{ w0  Wj  wq'  and  v^'J ,  the  curve  reduces  to  a  conic. 

We  have  the  equation 


wv  =  [  wx  wy  w  J 

=  fu^  u2  u  1]  A. 

In  this  case,  A  is  a  4  x  3  matrix.  Now  if  the  top  row  of  this  matrix  is 
( 0  0  0  J ,  the  equation  reduces  to 
2 

wv  =  [  u  u  1  ]  A 


when  the  top  row  of  A  has  been  omitted.  A  is  now  a  3  x  3  matrix,  and  it  is 
possible  to  show  that  this  equation  is  a  parametric  form  for  the  general  conic, 
expressed  as  a  quadratic  rational  function. 

For  the  top  row  of  A,  we  have  the  vector  equation 


{2-2  11] 


Expanding: 


wovo 

W1 V1 

V  vo 

+ 

woV 

V  V1 

+ 

W1  V 

[0  0  0) 


2  W0  vo  ”  2  W1  V1  +  V 


vo  +  wovo  +w 


V1  +  W1  V1 


t  _ 


=[000] 


Collecting: 


w0 <2v0 +  V)  +  W1 ("2 V1 +  Y> +  w0f  <V +  V <V1>  =  f°  0  °i 


In  matrix  form. 


[0  0  0]  =  |w0  Wj  w#'  Wj’]  r‘ 


2  vo  +  V 


-2  vx  +  v  ' 
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The  matrix  is  a  4  x  3. 

Now  we  can  adjoin  a  column  to  the  matrix  and  an  element  to  the  vector, 
and  write 


[0  0  0  w0I  «  fw0  wx  w0'  w^J 

2  v0  4  V 

1 

-2  Vj  4  », 

0 

vo 

0 

.  V1 

0 

then 


fw0  wi  wo'  W1*  j  *  (  0  0  0  w0j 

'2  v0  4  vo' 

1 1 

-2  v  +  v  ' 

A  * 

c 

vo 

0 

.  V1 

0 

m 

If  the  indicated  inverse  exists,  then  a  solution  can  be  obtained  in  terms  of  wQ. 
Furthermore,  wQ  can  be  set  equal  to  1  arbitrarily. 

The  matrix  has  an  inverse  in  case  the  determinant 


-2  Vj  +  v  * 


^  0. 


As  a  test,  construct  a  conic  with  end  conditions 


V 

'oof 

V1 

1  1 1 

V 

10  0 

.  V. 

10  0 

v,*- 


v 


1 


I 
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2  v0  +  V 

[0  0  2]  +  (1  001 

ss 

1  0  2~ 

-2  Vj  +  Vj' 

( -2-2-2  1  4(100] 

-1  -2  -2 

V0 

0  0  1 

0  0  1 

.  V1  . 

111 

111 

» 

We  first  test  to  see  whether  the  determinant  vanishes: 

1-1  -2  -2 

I  1-1  -2  I 

=  - 1.  Hence  the  augmented  matrix 


0  0  1 
111 
will  have  an  inverse. 


-1  -2 

1  1 


The  matrix  is 

10  2  1 

ami  its  inverse  is 

- 1 

0 

M 

O 

to  [ 

-1  -2  -2  0 

0  -1  -1  -1 

0  0  10 

0  0  10 

O 

rH 

fH 

*•4 
_ 1 

1  -1  -2  -2 

then 


[1  Wj  w0'  Wi'l  =  [°  0  0  i) 


0  10  2 

0  -1  -1  -1 
0  0  10 

1  -1  -2  -2 


whence  {1  w^’  ]  = 

The  conic  equation  ie 

wv  «  { u3  u2  u  1  ] 


(1-1-2  -2], 


2-211 
-3  3-2-1 
0  0  10 
10  0  0 


10  0  0 
0-100 
-2010 
0  -2  -1 


0  0  1 
111 
10  0 
10  0 


■-  A- 
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=  [  U3  U2  U  1 1 


1-1-2  1 
-2010 
10  0  0 


0  0  0 

-2  -1  0 

10-2 


■ 

r~t 

o 

o 

L 

111 

10  0 

J 

_1  0  0_ 

-  1 

u2  u  1] 

-2  -1  0 
10-2 
0  0  1 


0  0  1 

[wx  wy  w]  »  [-2i*2  +  u  j-u2]  -2u  +  l] 


X  =  u 


2 

-u 

y  -  *T^2u" 

asymptote  at  1  -  2x  =  0, 


x  = 


The  curve  is  hyperbolic,  with  an 


It  is  always  possible  in  all  of  the  foregoing  to  set  ~  1.  This  is 
because  all  equations  are  homogeneous.  It  is  never  possible  for  \vQ  =  0, 
since  this  leads  to  certain  degenerate  cases. 


We  remark  in  passing  that  when 


(W0  W1  V  V1  =  !1  1  0  °1 

the  equation  reduces  to  the  ordinary  parametric  cubic,  given  by 


v  =  [uVul]  M 


and  w  is  constant  and  equal  to  1.  Hence  the  rational  polynomial  functions 
contain  as  special  cases  all  conics,  ordinary  cubics,  and  of  course  therefore 
straight  lines  and  circles. 


* 
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Their  use  as  boundary  curves  for  surface  [>atehes  is  obvious.  They 
maintain  tangent  vector  continuity  between  adjacent  patches;  indeed,  if  the 
Fq  and  functions  are  constructed  as  rational  functions,  we  can  establish 
the  Fj  function 


^(u) 


3u  -  3u  +  1 


This  function  has  the  end  conditions 


V 

0 

r- 

0 

0 

1  ' 

V1 

1 

1 

1 

V 

= 

1 

0 

0 

V 

1 

0 

0 

V 

0 

0 

0 

.V. 

0 

0 

0 

Since  v^"  =  Xv0'  =  0(1  0  0]  and  Vj"  =  Xv^  =  0(1  0  Oj,  the 
curve  has  a  point  of  inflection  at  u  =  C  and  u  =  1.  Hence  its  use  insures 
curvature  continuity  across  boundaries  between  patches,  provided  of  course 
the  boundary  curves  have  similar  curvature  continuity  at  patch  corners. 

The  cure  is  symmetric.  Furthermore,  we  can  put 

Fq(u)  =  1  -  F^u) 

and  obtain  directly  the  Fq  function,  another  cubic  rational  function,  with 
similar  properties  to  F^. 

8.3  AN  EXAMPLE 

We  shall  work  out  the  equation  for  the  F,(u)  blending  function  with  the 
customary  stipulations  that 

Fjfl)  =  1,  Fx(0)  =  F  1*(0)  =  Fj'fl)  =  0, 

and  with  the  two  additional  stipulations  that  Fj,"(0)  =  F^"(l)  =  0  as  well. 
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I 

I 

i 


i 

i 

i 


l 

i 


This  blending  function  will  give  both  slope  and  curvature  continuity  across  the 
common  boundary  between  two  contiguous  patches.  The  end-conditions  are, 
forv  -  [F0(u)  1], 


V 

"0 

1 

V1 

1 

1 

V 

0 

0 

V 

0 

0 

V 

0 

0 

V  ,f 

L  1  J 

0 

0 

The  matrix  [  P  Q  R  j 


V"  0  0  o" 

0  V  "  0  0 

0  0  0  0 

.0  0  0  0. 


0  -6  0  4~ 

6  0  4  4 

0-201 
-2  0  -1  -1 


obtained  by  direct  substitution  in  the  given  form. 


Us  inverse  is 


0-1  0-4 

10-40 
-2  2  6  6 
2  0-60 


=  S 


5 


[ 


) 

1 
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Now  set  v„  **  — 

w  -  Z 


[4-  xl 


for  svmmBtrv.  Then  since  a  = 

V  - «#  •  * 


[W0  wx  w0'  Wj'l  =  [0  0  8  vcl  S 


[0  0 


4  8]  I" 


0-1  0-4 

1  0  -4  0 

-2  2  6  6 

2  0-60 


=  (8  8  -24  24]. 


Now  if  Wg  a  1,  instead  of  8,  the  equation  becomes 


Jw0  W1  W0*  wi'!  a  I1  1  '3  3] 


we  have 


“  K 

w_ 


V  vo  4  wo  V 

W^  Vj  4  W^  V  * 


and  substituting  the  values  of  ( wQ  w,  w^’  w^1  ] : 


A  a  1  2  -2 


2  1  l]  To  1 

3-2-1  11 


-3  3 


9  0  10 


1  0  0  0  3  3 


0  -3 


1  0 
0  3 

0  -3 


0  1 


Finally,  wv  =  |wxw]  =  [uu“ul]  1  0 

0  3 

0  -3 
0  1 


r  wx  w  ]  =  [  u3  |  3u2  -  3u  4  1  ] 
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Hence  x  = 


Fg  function  is 


w  3u2  -  3u  +  1 


=  F^(u)  as  required.  The  other 


F0(u) 


1-F  (u)  =  1 - - - 

3u  -3u  +  1 


3  2 

-u  +  3U  -  3u  +  1  . 

3u2  -  3u  +  1 


This  Is  seen  to  be  a  rational  cubic  function  also. 

8. 4  PLANE  CURVE  THROUGH  A  POINT 

The  plane  curve  vector  is  v  f  x  y  I  j ,  wv  -  [  wx  wy  w  ].  Here 

wx  wy 

the  polynomial  denominators  i  4  x  =  —  and  y  =  — ■*-  are  both  the  same. 
J  w  w 

As  before,  the  end  conditions  on  the  curve  are  contained  in  the 

matrix  fv„  1  . 


We  wish  to  cause  the  curve  to  pass  through  some  arbitrary  point  v 

C 

(commonly  called  a  "shoulder  point")  and  it  will  turn  out  that  we  shall  also 
be  free  to  choose  some  arbitrary  slope  at  this  point.  It  is  important  to 
distinguish  between  the  term  "slope"  and  "tangent  vector".  The  slope  of  a 
curve  implies  that  the  direction  of  the  tangent  vector  is  known,  but  the  magnitude 
of  the  vector  is  not  under  our  control. 


We  begin  by  assuming  some  value  of  the  parameter  u  to  correspond  to 

v  .  For  purposes  of  illustration,  let  u  =  at  this  point. 

C  4 
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i hen,  from  wv  =  [u'u"ul)  A,  we  write 


v  = 
c 


fl  2  4  8]  M  w0vQ 


W1V1 

V  V0  +  wo  V 


w  '  V1  +  W1  V1 


Observe  that  we  have  arbitrarily  set  w  =  1.  This  is  harmless,  since 
the  equation  is,  as  we  have  observed,  homogeneous. 

By  multiplying  the  matrices,  combining,  and  collecting  terms  as  we 
have  done  before,  we  achieve  the  result 

8VC  =  (W0  Wj  w0'  w  >1 


4  v  -  v  1 
1  1 


The  matrix  on  the  right  is  a  4  x  3  matrix;  to  make  it  square,  so  that  it 
can  have  an  inverse,  we  need  an  additional  column.  This  column  can  be  pro¬ 
vided  by  a  scaler  equation,  and  the  slope  relationship  will  furnish  this  equa- 


We  first  find  an  expression  for  the  tangent  vector  at  vc.  Differentiating, 
we  obtain  as  usual. 


(wv)'  =  [  3u  2u  1  0  J  A 

w*  *  [  3u2  2u  1  0  ]  M  f  w. 
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This  last  equation  comes  from  the  equation  for  (wv)'  by  an  argument  that 
we  have  used  before;  that  since  wv  =  [  wx  \vy  w  j ,  (wv)'  -  [  (wx)'  (wy)'  v/' ) 

Hence  w'  corresponds  to  the  last  column  of  fw  v 

1  o  o 

i  W1  V1 


<wovo>1 


which  is  f  w 


Now  at  v  ,  we  have  already  set  u  = 


We  make  this  substitution,  and 


obtain,  from 


wv’  =  (wv)'  -  w'v 

v  **r  [■*  *  -i  *i  \  <\  -  v 

*1  (v,  -  *c> 

"o'  (vo  -  V  +  wo  vo'  • 

V  <vi  -  v  +  »i  V 

When  we  perform  the  indicated  multiplications,  and  then  collect  results 
and  restore  to  matrix  form,  we  have 


4  v  '  = 
c 


'  -  i»0  vi  V  V  [-« (v0  -  V  -  V 


6{,1'VC>-V 


v  -  v„ 
c  0 


v  -  v, 
c  1 
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We  are  now  ready  to  introduce  the  slope  condition.  We  could  write 

~  *x*  '  would  yield  awkward  results  when  the  slope  became  very 

great  and  approached  an  infinite  slope.  Instead,  we  choose  two  numbers  a 
and  b  so  that  ax'  »  by*.  It  is  obvious  that  these  numbers  can  very  appropri¬ 
ately  be  a  a  sin  6  ,  b  *=  cos  Q  where  Q  is  the  slope  angle.  Then,  for  v  , 

C 

•  t 

0  =  by  -  ax 
c  c 


This  is  a  scalar  equation. 

Now  x*c  is  the  first  component  of  v‘c,  and  y'c  is  the  second  component 

of  v'c.  Hence  these  quantities  correspond  tc  the  first  and  second  columns  of 

Hie  matrix  in  the  equation  for  v'  ,  respectively.  We  can  write  this  out  in 

c 

detail: 


0  = 


-6<y0-ycHr0' 


-a 


W "  V 


fWr 


Wj'l 


where  the  brackets  on  the  right  enclose  the  resulting  4x1  matrix  (or  column 
vector). 

The  factor  4  (of  4  v  ')  obviously  drops  out  of  the  equation, 
c 
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We  now  adjoin  this  to  the  equation  for  v  .  Again  we  can  drop  the  8  (from 

© 

8  v  ),  and  obtain 
© 


w. 


Wj’] 


4  V-  ♦  v  • 

0  0 

P 

1 

< 

*-* 

q 

V- 

r 

0 

-v. 

s 

L  1 

- 

The  vector  on  the  left  consists  of  four  components:  { xc  yc  1  0  J ;  the 
matrix  on  the  right  is  a  4  x  4  matrix.  If  it  has  an  inverse,  S,  then  we  can 
solve  for  [  w^  w^  w^'  w  '  J  by  the  equation 

rw0  w1  w0‘  w^J  =  {vc  0]  s. 


These  values  of  the  w  vector  cause  the  curve  to  satisfy  the  desired 
conditions. 


8.5  SECOND  DERIVATIVE  VECTORS 

We  have  already  discussed  rational  functions  for 
v  =  f  x  1 J  (and  of  course  f  y  1 J  and  f  z  1  ] . ) 


In  particular,  we  showed  that  curves  based  upon  these  functions  can 
usually  have  arbitrary  first  and  second  derivative  vectors  at  the  end-points, 
and  in  addition  can  be  caused  to  pass  through  some  shoulder  point  v^,  also 
arbitrarily  chosen. 

When  the  vector  v  =  [  x  y  1  ] ,  the  complete  generality  of  the  resulting 
curves  is  somewhat  curtailed.  We  shall  investigate  the  conditions  under  which 
such  a  plane  curve  can  satisfy  end  conditions  including  second  derivatives. 

We  have  already  obtained  an  expression  for  the  vector 

(woV  *ivi"l  =  two  wi  w0’  wi*  1 

Before,  the  vectors  wQ  vQ"  and  w^  v  ”  were  actually  scalars,  since  they  came 
from  v”  =  { x”  0  J .  But  since  v"  =  ( x”  y’!  0  J ,  they  are  each  2-component 
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vectors,  and  their  combination  makes  a  4-component  vector.  Similarly  P  and 
Q  are  now  each  a  4  x  2  matrix,  and  their  combination  is  a  4  x  4  matrix. 

We  carry  wQ  vQ”  and  v  "  across  the  equal  sign,  and  obtain  on  the  left 
the  null  vector: 

[0  C  0  01  «  [v’0  w2  w0'  Wj’J  /[P|Q]  "  "  v0"  0  0 

0  0  v  " 

0  0  0  0 

\  .0  0  0  0 

The  matrix  on  the  right  is  4  x  4.  Now  the  condition  that  must  hold,  in 
order  for  there  to  be  a  solution  for  [  wQ  wQ'  w  '  J ,  is  that  this  matrix 
must  be  singular;  the  determinant  of  this  matrix  must  vanish. 

This  last  remark  tells  us  that  v^"  and  v^"  cannot  be  chosen  entirely 
arbitrarily.  However,  it  is  always  possible  to  make  the  determinant  of  the 
matrix  vanish  by  the  adjustment  of  any  one  of  the  four  components  of  v^"  and 
v^".  Thus  if  one  of  the  four  components  is  the  number  a,  we  can  expand  the 
determinant  in  such  a  way  as  to  obtain  the  equation. 

k  a  +  k  -  0,  from  which  a  can  be  found. 

Suppose  the  matrix  is,  or  has  been  caused  to  be,  singular.  Then,  if 


[0  0  0  01  =  (wQ  wx  w w1'l  S, 

we  make  it  non-singular  by  an  appropriate  modification.  In  some  cases,  this 
might  consist  in  adding  1  to  an  element  in  the  top  row  of  S.  As  an  illustration, 
we  might  have 


[0  0  0  w0l  =  [w0  wQ’  w^l  S  + 


0  0  0  1 
0  0  0  0 
0  0  0  0 
0  0  0  0 
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ft  can  be  seen  that  this  modification  is  still  a  valid  equation;  if  the 
modified  S  matrix  now  has  an  inverse,  we  can  immediately  obtain  it  and  solve 
for  the  w  vector. 

As  an  illustration,  consider  the  end  conditions 


* 

-  — 

vo 

0  0  1 

V1 

111 

•» 

© 

► 

10  0 

V 

10  0 

V 

10  0 

vl" 

a  o  o 

L  X  J 

. 

We  plan  to  adjust  a  in  v  "  until  the  matrix  is  singular. 

We  require  first  the  matrix  [  P  |  Q  ]  ,  which  is  given  by 


—4  v  * 

4  vo 

«  <v0  -  vt)  .  2  v0’ 

6  (V,  -  v0)  -  2  Y 

4  v  * 

1 

-2  v  • 
c 

2  <v0  -  V 

-2  <v2  -  vQ) 

'2  V/ 

Substitution  of  the  end  conditions  gives  the  4x4  matrix; 


-4  0 

[-6  -6]  +  [2  0] 

-4 

0 

-4  -«1 

[6  6]  -  [2  0] 

4  0 

4 

6 

4 

0 

-2  0 

-2  -2 

-2 

0 

-2  -2 

-2  -2 

-2  0 

-2 

-2 

-2 

1 

°J 

In  passing,  we  note  that  if  our  end  conditions  had  been  v  n  =  [0  9  0] 
and  if  also  v  "  =  [0  0  C  ] ,  the  resulting  matrix  would  be  singular,  because 
the  first  and  third  columns  of  [  P  >  Q  ]  are  identical.  However,  this  is  not 
our  present  case. 
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We  now  subtract  the  matrix 


r 

«* 

m 

V 

0  0 

1  0 

0  0 

o 

o 

V 

0  0 

a  0 

0  0 

0  0 

0  0 

0  0 

.0  0 

0  0 

0  0 

0  0_ 

from  (  P  J  Q  ]  and  obtain 

'-5  0  -4  -6 " 

4  6  (4-a)  0 

-2  0  -2  -2 

-2  -2  -2  0 


By  a  series  of  reductions  accomplished  by  multiplying  rows  of  the  matrix 
and  additions  (or  subtractions)  of  rows  to  remove  elements,  we  can  obtain  the 
determinant 


<a-l)  -  1 
1  -1 


=  0 


This  implies 


a  =  2. 


This  is  the  value  of  a  that  makes  the  matrix  singular. 


The  singular  matrix  is 


-2 


0-4-6 
6  2  0 
0  -2  -2 


-2  -2  -2  0 


We  make  it  non-singular  by  adding  1  to  the  top  left  element  The  matrix 
is  now 

0  -4  4 
4  6  2  0 

-2  0  -2  -2 
-2  -2  -2  0 
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and  corresponds  to  the  vector  [  0  0  0  ]  on  the  left  of  the  equation. 

The  inverse  of  this  matrix  is 


~  2  1-63 

-10  3-1  =  R. 

-1  -1  3  -3 

-1  0  2  0, 

Finally,  [  w^  Wj  w  ’  w  *  ]  =  fwfi  0  0  0]  R. 

If  we  arbitrarily  set  wQ  =  2,  then  the  required  solution  is  just  the  top 
row  of  R,  or 

[W0  Wi  wQ*  w1'J  =  [2  1  -6  3}. 

With  these  numbers  known,  the  curve  equation  is  completely  defined.  We 
obtain  it  by  substitioo  in  the  canonical  form: 


W„  V- 

0  0 

W,  V, 

1  1 

w  *  v„ 

+  W_  V  » 

0  0 

0  0 

W  ’  V, 

+  W,  V,  ’ 

1  1 

1  1 

=  [u^^ul]  M  j~  0  0  2 


111 

[0  0  -6]  +  [2  0  0] 
f  3  3  3J  +  [1  0  0J 


3  2  .  . 

=  [  u  u  u  1  J 


2-211 
-3  3  -2  -1 

0  0  10 

10  0  0 


0  0  2 

111 
2  0-6 
4  3  3 
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=  [  U  1  1  4  1  -1 1 

-5  0  6 

2  0-6 

0  0  2 

The  separate  equations  for  the  x  and  y  coordinates  of  the  curve  can  he 
exhibited: 

«  3  _  2  _ 

4u  -  5u  +  2u 

X  ~  g  2 

-u  +  6u  -  6u  +  2 

3 

u 

y  =  — - - -  • 

-u  +  6u  -  6u  +  2 


In  the  foregoing,  certain  matrices  have  occurred.  These  matrices  are 
significant  ones,  and  can  be  written  as  transformations  of  the  common 


as  follows: 


For  the  conic  condition  matrix. 


t 


! 

1 
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For  the  P  and  Q  matrices  associated  with  w  v  " 

0  0 


i 

i 


I 

I 

i 


! 
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APPENDIX 

COMPUTER -GENERATED  IMAGES 


Three  orthographic  views  and  a  perspective  view  (upper  right)  of  a  single 
surface.  The  perspective  was  generated  by  the  computer  from  the  three  orthog¬ 
raphic  views  given  it.  When  registered  on  a  cathode-ray  tube  a  change  in  any 
one  of  the  views  will  automatically  cause  a  change  to  be  made  in  the  other  views. 

•These  computer-generated  images  were  provided  by  Prof.  B.  Herzog  of  the 
University  of  Michigan,  and  are  shown  through  the  courtesy  of  Ford  Motor 
Company. 
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COMPUTER-GENERATED  IMAGES 


Perspective  views  of  other  surfaces. 
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